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PURPOSE : To analyze the long-term behavior of earth satellites due to the gravita- 

tional effects of the sun, the moon, and the oblate earth. 

RELATED TO : RAND's continuing study of satellite orbit control requirements for NASA, 

The study extends the results obtained in R-399-NASA, Perturbations of a Synchro- 
nous Satellite , May 1962. 

DISCUSSION AND METHODOLOGY : Now that long-lived earth satellites are orbited on a 

relatively permanent basis, it is important to know the effect of long-term or- 
bital perturbations, and the cost of controlling them. It is also interesting 
to consider the effect of such perturbations on the motion of the increasing 
debris in orbit. While the earth's inverse-square law gravitational field is 
the major attraction on a satellite, perturbing forces arise from the attraction 
of the sun, the moon, and the oblateness of the earth. The resultant force field 
is neither central in direction nor inverse-square in magnitude. The components 
of the perturbing forces that are normal to the orbital plane perturb the plane's 

orientation relative to inertial space. Those 
lying in the orbital plane cause changes in 
the shape and orientation of the orbit in its 
plane. The in-plane motion was analyzed in 
R-399— NASA; in the present study en^hasls is 
on the determination of the orbital plane it- 
self. The analysis applies to satellites in 
near-circular orbits at any inclination and 
with orbital radii less than 10 earth radii. 

The perturbed motion of an uncontrolled satel- 
lite is described as seen from inertial space 
and as seen from the rotating earth. 

PRINCIPAL FINDINGS : The motion of the orbital 

plane is such that its normal describes a coni- 
cal surface relative to inertial space as shown 
in the figure. The ground trace of a synchro- 
nous-altitude orbit lying in the reference 
plane is a figure eight with crossing-point on 
the equator and a maximum latitude excursion of 
7 20’; this does not vary with time. An orbit 
at an angle to the reference plane has a figure- 
eight ground trace which varies with the re- 
gression period. For a synchronous orbit that 
is originally equatorial and "stationary," the 
ground trace develops from a point to a figure 
eight, with a latitude excursion of 14°40' after 26,6 years, and then reverses the 
process. The total regression period is about 53 years, A fuel expenditure pro- 
portional to the sine of twice the inclination angle relative to the reference 
plane is required to stop orbital regression. 
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PREFACE 


As part of a continuing RAND study for the National Aeronautics 
and Space Administration of the perturbations of earth satellites and 
resulting orbital control requirements, this Report investigates the 
gravitational effects of the sun, the moon and the oblate earth on the 
orbital behavior of earth satellites. An extension of the results 
obtained earlier in R-399-NASA, the analysis provides determination of 
orbital control propulsion requirements and of the orbital regression 
of earth satellites with no restriction on orbital altitude or incli- 
nation. The general analytical solution for the regression of a satel- 
lite obtained in the R can also be specialized to explain lunar re- 
gression and the behavior of low-altitude satellites. 

The Report should be of interest to people dealing with satellite 
systems which require high-precision long-term orbital control. 





SUMMARY 


This Report presents an analysis of the long-term behavior of earth 
satellites due to the gravitational effects of the sun, the moon and the 
oblate earth. The importance of these effects has increased in recent 
years because of interest in more precise orbital control and an increase 
in expected payload lifetime. The fuel requirement for maintaining pre- 
cise long-term orbital control in the presence of these gravitational 
perturbations tends to dominate the overall fuel requirement for orbit 
and attitude control. As a result, it is of increasing interest to de- 
termine first whether the magnitude of the long-term orbital perturba- 
tions of an uncontrolled satellite is compatible with the requirements 
of a given satellite mission during its expected payload lifetime. If 
the orbital perturbations exceed acceptable values , it is necessary to 
investigate methods whereby these perturbations can be controlled, 
either actively or passively. In addition, it is of interest to con- 
sider the effect of these long-term perturbations on the motion of the 
increasing collection of debris in orbit. This Report deals with all 
of these problems . 

The analysis presented here applies to satellites in near-circular 
orbits at any inclination and with orbital radii which are small compared 
to the radius of the moon's orbit (i.e., less than 10 earth radii). The 
formulation takes into account the rotation of the earth around the sun, 
the rotation of the earth-moon system about its center of mass and the 
regression of the moon's orbital plane about the normal to the ecliptic. 

The major effect of the perturbing influences considered is to pro- 
duce motion of the orbital plane relative to inertial space. The nature 
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of this motion can be completely described by the trace of the normal 
to the orbital plane on a sphere concentric with the earth. It is shown 
that for an orbit of a given radius an orbital orientation can be found 
which remains invariant relative to inertial space. This invariant plane 
has a common intersection with the earth’s equatorial plane and the plane 
of the ecliptic, while its inclination to the latter is always less than 
that of the equatorial plane. For low-altitude orbits, the invariant 
plane is very nearly equatorial, with an inclination of 23°27' relative 
to the ecliptic. As the orbital altitude increases, the value of the 
inclination decreases to 16°7' at synchronous altitude and approaches 
zero for extremely high orbits. 

For an orbit of a given altitude, the trace of the normal to the 
invariant plane on the sphere described above is a single point between 
the earth's polar axis and the normal to the ecliptic. For orbits of 
the same radius but different orientations, two types of motion are pos- 
sible. If the initial inclination of the orbit relative to the corre- 
sponding invariant or reference plane is less than about 80^, the normal 
to the orbital plane rotates about the normal to the reference plane with 
an essentially constant angular rate and inclination angle. The resulting 
trace on the sphere is a circle with center on the normal to the reference 
plane. If the initial inclination is in excess of 80°, the trace of the 
normal to the orbital plane on the sphere may be an ellipse with its 
center on a line in the direction of the vernal equinox and major axis 
in the reference plane. 

The regression period at zero inclination varies from .1 year for 
a surface orbit to about 53 years for a synchronous orbit and a maximum 
of about 75 years for an orbit of radius equal to 9 earth radii. As the 
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inclination increases, the period varies inversely as the cosine of the 
inclination angle. On the other hand, the regression period corresponding 
to the elliptical contours has a minimum value as the ellipse approaches 
a point, and increases toward infinity as the major axis approaches 90°. 
However, since high-inclination orbits are of relatively little interest, 
the emphasis in this Report is on the first type of regression, which is 
illustrated by the frontispiece. 

It should be noted that superposed on this steady-state motion are 
oscillatory perturbations in both regression rate and orbital inclination 
which cause the instantaneous position of the normal to the orbital plane 
to oscillate relative to its steady-state motion. However, it is shown 
that the displacement is less than half a degree. 

It is of particular interest to observe the effect of this orbital 
regression on the relative motion of synchronous altitude satellites as 
seen from the rotating earth. Since the orbital altitude is assumed to 
be constant, this relative motion is completely described by the trace 
of the subsatellite point on the earth's surface. 

The reference plane corresponding to a synchronous altitude orbit 

o 

has an inclination of 16 7 relative to the ecliptic, as compared with 

an inclination of 23°27' for the earth's equatorial plane. Since the 

orientation of a synchronous orbit in this plane remains invariant rela— - 

o 

tive to inertial space, its inclination of 7 20 relative to the earth s 
equatorial plane is also invariant. As a result, the trace of the sub- 
satellite point on the surface of the rotating earth is the character- 
istic figure-eight pattern with a maximum latitude excursion on either 
side of the equator equal to the inclination angle of 7°20' relative to 
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the equatorial plane. In addition, the maximum longitude excursion 
relative to the equatorial crossing position is of the order of ±15 min 
of arc. This ground trace is repeated once each orbit with no variation 
in size or shape. Similarly, if the orbital plane of a synchronous orbit 
is perpendicular to the reference plane and polar relative to the earth, 
it remains stationary relative to inertial space, and its ground trace 
on the rotating earth also repeats itself on each orbit. However, for 
such an orbit, the ground trace varies from -90° to ±90° in latitude 
each day. 

For any other inclination of a synchronous orbit relative to the 
reference plane, it is found that the maximum latitude or amplitude of 
the figure eight varies as a function of time. This is due to the fact 
that the maximum latitude is equal to the inclination of the orbit rela- 
tive to the equatorial plane, and although the inclination relative to 
the reference plane is fixed, that relative to the equatorial plane varies 
as the orbit regresses. The resulting variation in the ground trace am- 
plitude has a periodicity equal to that of the regression and a magnitude 
which can never exceed 14°40' . In addition, it is found that the longi- 
tude of the equatorial crossing also oscillates with the regression 
period and with an amplitude which may be as large as 7°20' , depending 
on trie orbital inclination relative to the reference plane. 

In regard to these variations in the size and shape of the ground 
trace, it is of particular interest to consider the long-term behavior 
of a satellite which is initially in a synchronous equatorial orbit. 

Such a satellite is ordinarily referred to as stationary since it appears 
to be fixed relative to the earth. However, its orbit is actually in- 
clined to its reference plane at an angle of 7°20* and has a regression 
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period of about 53 years. As the orbital plane regresses, its inclina- 
tion to the equatorial plane increases from 0*^ at an average initial 
rate of .863° per year. At the end of half the regression period this 
inclination reaches a maximum of 14°40', after which it decreases sym- 
metrically to 0° after a complete regression period. Since the maximum 
latitude excursion during an orbit is equal to the orbital inclination 
to the equator, the resultant ground trace is initially an equatorial 
point but develops into a figure eight which reaches a maximum amplitude 
of 14°40' before decreasing to the original equatorial point at the end 
of the regression period. During this cycle, the position of the equa- 
torial crossing oscillates with the regression period with an amplitude 
of .47°, moving to the east of its initial position during the first 
half of the cycle and to the west during the second half. It should be 
noted that these variations in longitude are superposed on the shorter 
period (~ 2-year) oscillations due to the equatorial ellipticity de- 
scribed in Ref. 1. 

It is seen that a passive satellite cannot remain truly stationary 
relative to the rotating earth, and that its earth trace can remain in- 
variant only for certain orbital inclinations. Since a given satellite 
mission may require a fixed ground trace which is not inherently in- 
variant, it is of interest to determine the amount of control necessary-— 
to produce the desired invariance. It is seen that such an invariant 
ground trace can exist only if the orientation of the orbital plane re- 
mains fixed in inertial space. By applying appropriate control impulses 
normal to the orbital plane, it is possible to reduce the steady-state 
orbital regression rate to zero. Under these conditions, the orbital 
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orientation and the resulting ground trace have the desired invariance. 
The control impulse required per year to achieve this invariance is pro 
portional to the sine of twice the inclination angle of the orbit to it 
reference plane. Thus, the magnitude of the control impulse per year 
depends as follows on the desired value of the ground trace amplitude. 
For an amplitude of 0°, the required control impulse has an average 
value of 152 ft/sec/year, which decreases to zero for an amplitude of 
7^20', after which it increases to a maximum of 580 ft/sec/year for an 
amplitude of 45°. The impulse requirement for amplitudes between 45° 
and 90° is a mirror image of that from 0° to 45°, decreasing to zero 
at 82°40' and increasing again to 152 ft/sec/year at 90°. It should 
be noted that these values represent a long-term average control re- 
quirement, neglecting the oscillatory components of the orbital re- 
gression. In the event that it is necessary to control these oscilla- 
tory variations, the control requirement in a given year might deviate 
from its average value by as much as 30 ft/sec, depending on the ampli- 
tude and phase of the oscillatory terms. 

If instead of an invariant ground trace, a given mission requires 
merely an upper limit on its latitude excursion, it may be possible to 
satisfy this condition passively. If the initial orbital inclination 
to the equator is made equal to the upper limit of the latitude excur- 
sion, then by a suitable choice of the initial regression phase the in- 
clination to the orbital plane will decrease to 0° before it again in- 
creases to its initial value . 

In this way the time during which the latitude excursion remains 
below its upper limit is maximized. If this time is longer than the 
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expected payload lifetime, this passive technique can be used to satisfy 
the mission requirement. However, the decision regarding the use of 
active or passive orbit control depends on the tolerances in ground 
trace amplitude and the required mission lifetime. 

Finally, the analysis of orbital regression developed here for 
artificial satellites is extended to include the regression of the moon. 
This requires an expansion of the basic theory to include higher order 
terms as shown in Appendix E . 
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I . INTRODUCTION 


In the ten years since the launching of the first Sputnik, a rather 
impressive amount of hardware has been placed in orbit around the earth 
in the fulfillment of various satellite missions. During the early years 
of this period, the lifetime of a given satellite mission was short and 
there was little interest in precise long-term orbital control. However, 
in recent years, as technology has developed which can take advantage 
of precisely controlled synchronous orbits with long-life potential, it 
has become important to know (1) the effect of long-term orbital per- 
turbations on the satellite mission, and (2) the cost of controlling 
these perturbations . 

It is well known that if the only force on a satellite is the in- 
verse-square-law gravitational attraction of the earth, then the re- 
sulting orbit is an ellipse with one focus at the earth's center. In 
addition, the direction of the normal to this orbital plane remains 
fixed relative to inertial space. While the earth’s inverse -square field 
is the major attraction on the satellite, it is necessary to consider 
the perturbing forces which might produce long-term changes in the basic 
orbital motion described above. Three such forces are those due to the 
attractions of the sun and the moon and that arising from the oblateness 
of the earth. Since the resultant force field when these effects are in- 
cluded is neither central in direction nor inverse -square in magnitude, 
the resulting orbital perturbations may be of two types. Those components 
of the perturbing forces which are normal to the orbital plane produce 
perturbations in the plane's orientation relative to inertial space. 

Those components which lie in the plane cause alterations in the shape 
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and orientation of the orbit in its plane. 

In Ref. 1, the effects of the sun and the moon on a satellite in 
a synchronous equatorial orbit are determined. The results show that 
the in-plane perturbations are of the nature of small amplitude oscil- 
lations in the satellite's position relative to its nominal unperturbed 

position. The maximum excursion is of the order of 45 mi. It is also 

found that the perturbations in the attitude of the orbital plane are 

of the nature of a slow change in its inclination to the equatorial 
plane at a rate of about .85 /year. The analysis also indicates a slow 
sinusoidal increase to a maximum inclination of about 20° and a return 
to 0 after a period of about 73 years. However, these two values are 
only approximate, since a 20° angle exceeds the small angle assumption 
used in the perturbation analysis. 

This Report gives a more general determination of the orbital per- 
turbations resulting from the gravitational attractions of the sun and 
moon and from the oblateness of the earth, and it places no restriction 
on the magnitude of the orbital inclination. The emphasis is primarily 
on the determination of the motion of the orbital plane since it is not 
anticipated that the in-plane motion will differ greatly from that de- 
termined in Ref. 1. By means of the analysis presented here, the per- 
turbed motion of an uncontrolled satellite is described both as seen from 
inertial space and as seen from the rotating earth. In addition, the fuel 
requirement to maintain a fixed orientation of the orbital plane relative 
to inertial space is determined. 
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II. METHOD OF ANALYSIS 
STATEMENT OF THE PROBLEM 

The problem to be solved in this Report can be stated as follows: 

If a satellite is in a circular orbit around the earth with a known 
initial orientation of its orbital plane relative to inertial space, 
determine the motion of this orbital plane as it is affected by the 
gravitational attraction of the sun, the moon and the oblate earth. 

For the purposes of this analysis it is assumed that the orbital 
altitude is sufficiently high that forces due to residual drag can be 
neglected. In addition, by assuming a small area-to-mass ratio for the 
orbiting object, the effects of solar radiation pressure can also be 
neglected. The ellipticity of the earth's equatorial section is also 
omitted since its effect on orbital regression is negligible. 

Finally, the positions of the sun and the moon relative to the earth 

are specified as known functions of time according to the following model. 

The center of mass of the earth-moon system moves on a circular orbit 

around the sun with a period of one sidereal year, and the plane of the 

"k 

motion is that of the ecliptic. The earth and moon rotate about their 
common center of mass with a constant separation and a period of one 
sidereal month. The plane of this rotation is inclined at an angle of 
5°8' to the ecliptic and regresses about the normal to the ecliptic 
with a period of 18.6 years. 


Strictly speaking, the plane of the ecliptic is defined by the 
motion of the earth's center of mass; however this differs in orienta- 
tion by about 1 sec of arc from the plane defined above. 



4 


DEFINITION OF COORDINATES 


Reference Systems 

In the formulation of the equations of motion and the description 
of the resulting motion, it is convenient to define the following ref- 
erence systems. 

X , y , z . This is an inertial reference system with its origin 

O Q i 

at the center of the sun; its x y plane is the ecliptic and its x axis 
is in the direction of the earth at the time of the vernal equinox. 

x^, y^, z^. This and all the remaining reference systems are geo- 
centric. This particular one maintains its axes parallel to the corre- 
sponding ones in the x , y , z system. Thus, every point of this sys- 

d O D 

tem is under a constant acceleration in a direction parallel to the earth 
sun line. The unit vectors associated with this and the previous sys- 
tem are represented as i , j , k . 

^ o o o 

x^ , y- , z, . This system is rotated relative to x , y , z through 
ill o o o 

an angle about their common x axes as shown in Fig. 1. This system 
is the one relative to which the motions of the orbital plane are ex- 
pressed. The appropriate value of is determined in the course of the 
analysis. The associated unit vectors are represented as i^, k^. 

X - , y - , 2 . This system has a common x axis with the two previous 

systems, while the x y plane is the earth's equatorial plane, which 

G G 

makes an angle X = 23^27' with the ecliptic or x^y^ plane as shown in 

Fig. 1. The associated unit vectors are represented as i , j , k . 

G G G 

X , y , z . This is the orbital coordinate system with the x axis 
passing through the satellite and the xy plane representing the instan- 
taneous orbital plane. The orientation of this system relative to the 




Fig.l — Reference systems 
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^1’ ^1’ specified by the three angles a, i|f and 9 as shown 

in Fig. 2. The angle o is the angle between the xy (orbital) plane 
and the x^^y^ (reference) plane. The angle \1; is the regression angle 
measured from the x^ axis to the line of nodes, ON. Finally the angle 
0 is the orbital angle measured from the line of nodes to the x axis . 
The associated unit vectors are represented as i, j, k. 

It should be noted that an additional relation between O', ijr and 
0 is required to insure that the xy plane is indeed the instantaneous 
orbital plane. This relation is determined in the derivation of the 
equations of motion. 

X , y , z . This is the lunar reference system in which the x 
m m m m 

axis passes through the center of mass of the moon and the plane 

represents the moon's orbital plane. The orientation of this system 

relative to the x , y , z system is specified by the three angles O' , 
o o o m 

\lf and 0 as shown in Fig. 3, where it is seen that these angles are 
m m 

analogous to o', ijf and 0 for the satellite orbital system. As before, 

the associated unit vectors are represented as i , i , k . 

m m m 

The direction cosines relating these various systems are listed 
in Appendix A. 

Position Vectors 

The relative positions of the sun, moon, earth and satellite can 
be described vectorially as shown in Fig. 4, in which the various vec- 


tors are defined as follows. 
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Fig. 2 — Satellite orbital reference systems 
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Fig. 3 — Lunar orbital reference system 
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The vector R is from the center of the sun to the center of mass 
of the earth-moon system. In accordance with the assumed motion of this 
system around the sun, 


R = 


Rr. 


R(i cos fi? + i sin 0) 
o o 


( 1 ) 


where R is constant and equal to the earth's orbital radius, while 0 
varies linearly with time. 

The vector ^ is from the center of the earth to the center of the 
o 

moon and lies along the x axis. It can be expressed as 

m 



( 2 ) 


where p is a constant and equal to the earth-moon separation, 
o 

The vector r is from the earth's center to the satellite along the 
X axis and can be expressed as 


(3) 


where r is the magnitude of r and is one of the orbital variables to be 
determined. 

In addition to these basic vectors, it is convenient to define cer- 
tain others that are of use in the formulation of the equations of 
motion, although they can be expressed in terms of those already de- 
fined. The vector p„ from the earth's center to the center of mass 

E 

of the earth-moon system can be expressed as 


I 

|J< 


P 


o 


(4) 
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where 


(1 = 




+ M 


m 


M 


( 5 ) 


m 


and Mg and are the masses of the earth and moon respectively. 

The vector from the sun's center to the earth's center is 
given by 


«E = ^ Po 


( 6 ) 


The vector R from the sun's center to the moon's center is given by 
m 


R = R + 
m 


(' - S Po 


( 7 ) 


The vector r from the sun's center to the satellite is expressed as 


r„=R-“p + r 

S M* o 


( 8 ) 


The vector r from the moon's center to the satellite is given by 
m 


r = r - p 
m o 


( 9 ) 


Angular Velocity Vectors 

The relative motion of the various coordinate systems defined pre- 
viously can be expressed in terms of the following angular velocity 
vectors . 

The angular velocity, , of the center of mass of the earth-moon 
system around the sun can be expressed as 
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<^ = 8 ( 10 ) 

where @ is a constant equal to .0172 rad/solar day. This corresponds 
to a period of one sidereal year. 

The angular velocity, o) , of the earth-moon system about its center 

m 

of mass and relative to inertial space is given by the expression 


au = 0 k + I k 
m mm mo 


( 11 ) 


where 0^ is the moon’s orbital rate of .22998 rad/solar day, corresponding 

to a period of one sidereal month, while t Is the regression rate of the 

m 

-4 

moon s orbital plane and has a value of -9.249 x 10 rad/solar day, 
corresponding to a period of 18.6 years. 

The angular velocity 0^ of the earth about its axis is given by 

\ = Vg 


where 0^ is equal to 6.3004 rad/solar day, corresponding to a period of 
one sidereal day. 

The orbital angular velocity, of the satellite around the earth's 

center relative to inertial space is expressed in terms of the three or- 
bital angles as follows 

% “ + Of (ij^ cos '1^ ^ 

where , a and 0 are the orbital variables to be determined together 


with r. 



13 


FORMULATION OF THE EQUATIONS OF MOTION 

The derivation of the equations of motion of the satellite relative 
to the earth is presented in this section in abbreviated form; the de- 
tails are given in Appendix B. 

The equations of motion of the satellite, the earth and the moon 
relative to inertial space can be expressed in vector form as follows 


Sate llite : 


r„ = R - — p 
S p. ^o 


+ r = 


F„ + F 
Ss ms 


+ F 


Es 


M 


( 14 ) 


Earth: 




R - - p 
p, o 


F + F 
SE mE 


«E 


( 15 ) 


Moon: 



( 16 ) 


In this derivation the dot notation for time derivatives signifies 
a derivative relative to inertial space, while d/dt is a derivative re- 
lative to the orbital (x, y, z) reference system. In addition, the sub- 
scripts S, E, m and s refer to the sun, the earth, the moon and the 
satellite, respectively. When they are used with M, the mass of the 
body they represent is indicated. When they are used with a vector 
force, F, the first subscript determines the attracting body, while the 
second defines the body on which the force acts. 
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The vector equation of motion of the satellite relative to the earth 
is obtained as the difference between Eqs. (14) and (15) in the form 


r = 


+ F + F^ 
Ss ms 

M 


F + F 
SE mE 


“e 


(17) 


An evaluation of the forces on the right reduces Eq. (17) to the follow- 
ing form 



.2 r 

m I — 3 N “ 

r - — (r . Pq) 

»o 


1 



(r • R) R 
R 


(18) 


where is the earth's equatorial radius, G is the universal gravita- ^ 
tional constant and is the coefficient of the oblateness term in the 
earth's potential function. 

The three component equations of motion corresponding to the x, y 
and z axes are obtained from Eq, (18) in the form 
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(r . 1) = - -3 [1 - 3 (kg . 1) J 

r 2r 


r0 


m r 


r — — 21 

1 - 3(i . ±y 

L m J 


- r(s)^[l - 3(r^ • i)^] 


(19) 


3J2GM^R; 


(r • j) (k^ • i) (k^ . j) 


3 r§‘ 


+ -TT^ O (i • j) 

M< m m 


+ 3 r 0 ^(r^ . i) (r^ • j) 


( 20 ) 


(r • k) 


3J_GM_R 

2 E o ,r Tn 


(k_ • i) (k^ . k) 


3r0 _ _ _ _ 

+ ' ^ (i • i) (i • k) 

|Ji m m 


+ 3re^(rj^ • i) (tj^ • k) 


( 21 ) 


The quantity r is given by the relation 


- d 
r 


~l + 2[u)^ + [“0 

at 
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from which the left sides of Eqs. (19), (20) and (21) can be expressed 
as 


(r • 


i) 






(23) 


(r • j) = 2 ^ . k) + r ■ k) + r(m^ ■ i) (u)^ . j) (24) 

(r • k) = - 2 ^ (cu^ . j) - ’ 3 J ^^o ‘ 


At this point, it is necessary to determine the constraint equation 
which insures that the xy plane is the instantaneous orbital plane. 
This condition requires that r as well as r must lie in the xy plane. 
This can be satisfied if 


(r 


k) = 0 


from which it follows that 


(u) • j) “ 0 

o 

and 



(26) 


(27) 


(28) 


Thus, Eqs. (23), (24) and (25) can be simplified to the form 


(r 


• i) 


d r 


dt 



• k) 


2 


(29) 
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(r * k) = * i) (^Q * ^) 

where 

(cjD^ • i) = \|f sin 9sino'+(icos 0 

((jj^ • j) = \^ cos Bsino'-orsin 0 

(uu^ • k) = ijf cos Q'+O 


(30) 

(31) 

(32) 

(33) 

(34) 


The complete equations of motion of the satellite are obtained 
by combining Eqs . (19) -(21), (27), (29) and (30) -(34) to give 


2 

- r(§ 4- \j; cos cv)^ = 
dt 


GMe 3J2®^'^o r - - 21 

- --^4 J 


m 

M- 


1 - 3(i 


m 




- r9^[l - 3(r^ • i)^] 


(35) 




3r0 _ _ _ _ 

+ — (i • i) (i • j) 


+ 3r@^(rj^ • i) (r^ • j) 


(36) 



18 


3J GM R _ _ _ _ 

a = - ^ 2 . 1 ) . k) cos 0 

(9 + ij/ cos Q')r 


30' 


+ 


m 


M<(0 + cos O') 


(i • i) (i • k) cos 0 


m m 


30' 


( 9 + ijr cos <y) 


(r^ • i) (r^ • k) cos 0 


(37) 


^ = - 


SJaGM^Ro 


(0 + tjf COS Q')r sin a 


(k^ • i) (k^ • k) sin 0 


3 0 


+ 


m 


|J,(9 + cos a) sin cy 


(i 


m 


i) (i • k) sin 0 
m 


+ 


30' 


(9 + ijf cos O') sin O' 


(r^ • i) (r^ . k) sin 9 


(38) 


These four equations constitute the desired orbital equations of motion 
in terms of the variables r, 0, or and ijf. 

MOTION OF THE ORBITAL PLANE — 

The exact solution of the equations of motion as presented above 
does not appear feasible since the terms resulting from earth oblate- 
ness, the sun and the moon are rather involved functions of the unknowns 
r, 9, (X and ijf. To simplify the problem, it is assumed that the quan- 
tities r, 0, Of and if can be represented as follows by a steady-state 
value plus a perturbation. 
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r = r + 6r (39) 

o 

9=0 +60 (40) 

o 

~ (41) 

a = a + 6(y (42) 

o 


Since the forcing terms due to earth oblateness, the sun and the moon 
in Eqs. (35) through (38) are in the nature of perturbations, it is 
assumed that in these terms r, 9, if and or can be replaced by their 
steady-state values, while 9 and ijf are expressed as 


0 = 0 t 
o 




(43) 

(44) 


In addition, it is assumed that 9 is much greater than if, and that 
the nominal orbit is circular. Thus, the steady-state values of r and 
9 are given by 

(45) 


r = r 


6 ^ - ^ 
r 

o 


( 46 >- 


With these assumptions, Eqs. (37) and (38) can be written as 
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2* 

3J R 9 _ _ _ _ 

. 1) . 1^) COS 9 


3 0' 


H (i- • i) (i • k) cos 0 

m m 


M.0 


30 ^ _ _ 

H ; — (r- • i) (r, • k) cos 9 

9 

o 


(47) 


2 • 

3J R 9 _ _ _ _ 

I — 2_-2— . i) gin 0 


r sin O' 
o o 


3 0' 


+ 


m 


M. 0 s in cy 
o o 


(i • i) (i • k) sin 0 
m m 


+ 


30' 


0 sin O' 
o o 


(r^ • i) (r^ • k) sin 0 


(48) 


thus giving or and ^ as functions of time. 

If the direction cosines tabulated in Appendix A are substituted 
in Eqs. (47) and (48), they can be reduced to the form 


127 

O' = ) A . sin (ju. t 

Z-i 1 1 


(49) 


i=l 


127 


t t r~ / B. cos u), t 

o s in cy Z_i 1 1 


(50) 


° i=l 
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The details of this transformation are presented in Appendix C together 
with a tabulation of the resulting expressions for the and 

values, as well as for ^ . 

An examination of this tabulation shows that the values are 

• • • • • 

linear combinations of the frequencies 0 , 0 , 0, i|r and ijr . In the 

o m m o 

case of A^ and B^, they may be functions of any or all of the quantities 

• * • 

a , a, , O' , 0 , r , J_ , © and 0 . The presence of J„ in a given ampli- 
tude indicates that this term is at least partially the result of earth 
oblateness. Similarly, © and 0 indicate contributions due to the sun 
and the moon, respectively. Thus, it is seen that the steady-state re- 
gression rate, as well as the first seven oscillatory terms are the 

result of a combination of all three perturbing influences. On the 
other hand, those terms for which 8 ^ i 22 are entirely due to solar 

influence, while the rest of the terms (23 i i ^ 127) are due to the 

moon. 

The solutions for or and f can now be obtained by integrating Eqs. 
(40) and (50) to give 

A. 

a = a + y “ r 1 - cos ou. t] (51) 

o Aj 00 ■ i 

i=l ^ 


127 


t t + 


sin O’ 


Y ^ 

Zj 00 , 


sin uo^t 


i=l 


(52) 


under the assumption that at t *= 0, or = and 0. 

It is now necessary to define the angle between the reference 
(Xfyf) plane and the plane of the ecliptic since the quantities if 
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and are all functions of (X^. The selection of is dictated by 
the fact that the analysis leading up to Eqs . (51) and (52) is based on 
the assumption that the angle ex between the orbital (xy) plane and the 
reference plane remains essentially constant and equal to O' . 

In view of this fact, a value of should be selected which mini- 
mizes the amplitudes of the oscillatory terms in Eq. (51). Particular 
emphasis should be placed on the low-frequency terms because of the 
factor 1/u)^. An examination of the A^ expressions in Appendix C shows 
that the only reasonable choice for is that value which makes A^ , 
the coefficient associated with the frequency, \|f^, identically zero, 
as follows. 

•2 2 

3 sin o^)J sin 2a^ - o^^) = 0 (53) 

@ r 

o 


1 + 


m 


2jj.0‘' 


(2 - 


which can be solved for in the form 


(54) 


This selection of a also makes the amplitudes A, , A^ , B , B and B 

■L 4 6 14 6 

identically zero since they all have the left side of Eq. (53) as a 
factor . 

In addition, the amplitudes A 2 > A^ , A^, B 2 , and B^ all involve 


tan 2(X, = 


•2 2 

2J„eV 

2 O Q 
2 

0 r 

o 


sin 2X 


* 2 •2 2 
e ^ 2 j.0^r^ 

1 ■* ^ (2-3 sin o ) + 2X 

2m. 0^ 0^r^ 


the factor f, given by 
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§2 

1 + — (2 ~ 3 sin^ O’ )1 sin^ or. 

r, m J 1 

2[jl@ 


•2 2 

2*^9 9 

, 2 o o . 2 ^ 

+ -.2 2 " 

0 r 


(55) 


While no choice 
termined by Eq. (54) 


of Q'^ can make the factor f vanish, the value de- 
does minimize f since 






(2 


3 sin' 



sin 2 <y^ 


- — 4°r sin 2 (X - a ) (56) 

© r 

o 

which is equal to zero from Eq. (53) . 

By the same type of reasoning it can be shown that this choice of 
or^ maximizes the amplitudes and associated with the frequency, 20^. 
However, this is a high-frequency term with a very small amplitude. 

Thus, it causes very little change in the inclination angle, ex. 

With the determination of it is now possible to use Eqs. (51) 
and (52) to describe the resulting motion of the orbital plane. In 
Appendix D it is shown the remaining oscillatory terms result in less 
than a .5° variation in the normal to the orbital plane. Thus, to a very 
good approximation, the motion of the orbital plane is represented by 
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or = O' 


t = tot 


(57) 

(58) 


which means that the orbital plane maintains a constant inclination, 
relative to the selected reference plane but regresses about the normal 
to the reference plane at a steady rate of given in Appendix C as 


3^^ cos a 


86 


fl + r (2 - 3 sin a )1 (2 - 3 sin or.) 

L o m J i 

2m.@ 


•2 2 

2J.0 R p « 

+ -W b - ^ - “i>] 


® r 


(59) 


Thus, the motion of the orbital plane is completely described by 
Eqs. (54), (57), (58) and (59). 

In obtaining the above result, it was assumed in Eq. (44) that t 
could be replaced by ^^t in Eqs. (37) and (38). However, an examination 
of Eqs. (50) and (59) shows that as the orbital inclination, cr^, approaches 
90°, the above assumption is no longer valid, since approaches zero 
and Eq. (50) is dominated by the oscillatory terms. ^ 

In Appendix F, the behavior of these high-inclination orbits is in- 
vestigated in more detail. It is found that the representation of the 
motion as described by Eqs. (57) through (59) is a good approximation of 
the actual motion as long as the regression takes place about the 


axis . 
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III. RESULTS AND DISCUSSION 

DETERMINATION OF ORBITAL REGRESSION 

On the basis of the analysis presented in the previous section, it 
is now possible to make a quantitative determination of the orbital re- 
gression rate and its axis of rotation. 

Reference Plane 

An examination of Eq. (54) shows that the inclination, , of the 

reference plane is a function of orbital radius. The relationship is 

shown in Fig. 5, which gives a as a function of r . It is seen that 

i o 

for low-altitude orbits, for which the earth’s oblateness is the dominant 
perturbing influence, the reference plane is very nearly coincident 
with the earth's equatorial plane = 23°27'). This is due to the 
fact that the equatorial plane is one of symmetry for the oblateness 
effect. Similarly, for high-altitude orbits, for which the combi- 
nation of the solar and lunar effects becomes dominant, the reference 
plane approaches the plane of the ecliptic = 0°). Again, this is 
due to the symmetry of these two effects relative to this plane. 

Steady -State Motion 

The steady- state regression of the orbital plane as specified by 
Eqs. (57) and (58) is represented in Fig. 6, in which the normal to the 
orbital plane (the z axis) maintains a fixed angle. O' , with respect to 
the normal to the reference plane (the axis). At the same time, the 
z axis rotates about z^ at an angular rate tracing a circular contour 

on the spherical surface as shown. The residual oscillatory terms 
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Orbital radius ( earth radii) 


Fig. 5 — Dependence of reference plane inclination on orbital radius 
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Fig. 6 — Steady-state orbital regression 
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in Eqs. (51) and (52) cause small variations in this trace, but, as 

shown in Appendix D, these residual oscillatory terms are negligible. 

While the z axis maintains a fixed angle, a , relative to the z 

o 1 

axis during its regression, it is seen from Fig. 6 that the angle a 

G 

between the z axis and the earth's polar axis (the z axis) varies be- 

G 

tween the limits 

Q-g = ho - - ap I ( 60 ) 

min 

and 


cy. 


max 


O' 

o 


+ (^ - 


(61) 


Thus, or^, which represents the inclination of the orbital plane to the 
earth's equatorial plane, will vary slowly as a result of the regression 
about the z^ axis. 

The regression period, the time required for the z axis to make one 
complete rotation about z^, is given by the relation 



(62) 


By combining Eqs. (54), (59) and (62), T can be determined as a function 

R ^ 

of the constant orbital inclination, , relative to the reference plane 
and the orbital radius, r^ . This relationship is shown in Fig. 7. It 
is seen that in general the period increases as the orbital inclination 
increases, becoming infinite for equal to 90°. The period also in- 
creases with orbital radius , but appears to reach a maximum in the vicinity 
of 10 earth radii. Presumably, this is due to the increase in the effect 



(SJD0X) pojjsd uoissajBa^ 



Dependence of regression period on orbital inclination and orbital radius 
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of the sun and moon and to the corresponding decrease in the effect of 
earth oblateness at this altitude. 

The numerical values of the regression period range from .099 year 
for a surface orbit at zero inclination to 52.84 years for a similar 
inclination at synchronous altitude. At a given altitude, the period 
for an inclination, , is obtained by dividing the zero inclination 
period by cos o'^ . 

Special Cases 

It is of interest to examine certain special cases of the expres- 
sion for the orbital regression rate given by Eq. (59). 

Earth Oblateness . If the oblateness of the earth is the only 

perturbing influence, its effect is given by setting S> and 0 equal 

m 

to zero in Eqs. (54) and (59). Under these conditions, Eq. (54) gives 
a value of equal to A.. Thus, the reference plane is the earth's 
equatorial plane, and the regression takes place about the earth's 
polar axis. Under these same conditions, Eq. (59) reduces to 

3J R^e 

^ ^ cos a (63) 

2r^ 

o 

which is the usual form for the regression due to oblateness. (See 
Ref. 2.) 

Sun and Moon . If the effect of the earth's oblateness is negligible, 
as in the case of large orbital radii, can be set equal to zero, with 
the result that the angle from Eq. (54) becomes zero. Thus, the refer- 
ence plane becomes the plane of the ecliptic, and the resulting regression 
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takes place about the normal to the ecliptic. The regression rate as 
given by Eq. (59) reduces to 



cos 


O' 

o 




3 sin 
8m- 0 

o 


O' ) cos O' 
m o 


(64) 


where the first term is the solar effect and the second is that due to 
the moon. Numerical evaluation shows that the regression rate due to 
the moon is approximately twice that due to the sun. 

Lunar Regression . The expressions developed in the previous sec- 
tion can also be used to determine the regression of the moon itself due 
to the influences of the earth's oblateness and the sun's gravitational 
attraction. As applied to the moon, Eq. (59) can be restated in the 
following form: 


30 cos O' 


m 


m 
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m 


•2 2 
J^e R 

0 p 

o 


- 3 sin' 


u] 


(65) 


where the first term in the bracket is due to the solar effect and the 
second term is that resulting from the earth's oblateness. Actually, 
the oblateness term is negligible, and Eq. (65) can be rewritten as 


30 cos O' 


m 


30' 


40 


m 


40 


m 


( 66 ) 


By means of this relation, the regression period of the moon is found 
to be 17.9 years instead of the accepted 18.6 years. The reason for 
this discrepancy is that Eq. (66) does not include the higher order 
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terms which are necessary in calculating the orbital regression rate of 
the moon. If these are included, Eq. (66) then has the form 



and the regression period based on this relation is the accepted 18.6 

years. The details of the derivation of Eq. (67) are given in Appendix E. 

While Eq. (67) is necessary in explaining the behavior of the moon, 

where the ratio @/0 is of the order of 1/13, the higher order terms are 

m 

negligible in the case of artificial satellites, for which this ratio 
is of the order of 1/365 or less. Thus, Eq. (66) is adequate for the 
purposes of this Report. 

APPLICATION TO SYNCHRONOUS ORBITS 

While the foregoing determination of orbital regression applies 
for any orbital altitude, the discussion will now be limited to syn- 
chronous altitude orbits . 

Condition for Synchronism 

From Eqs. (54) and (59), it is seen that for a synchronous alti- 
tude orbit (r^ = 26195.2 mi) the inclination of the reference plane, ^ 
is equal to 16°7' , while the steady-state regression rate is given 
by 


= 1(0) cos Oi 
o ^ ' o 


( 68 ) 
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-4 

where t (0) has a value of -3.257 X 10 rad/solar day. In the case 
of an unperturbed equatorial satellite, synchronism is achieved by set- 
ting the orbital angular rate, 0^, equal to the earth's spin rate relative 
to inertial space, 0 . Under these conditions the satellite remains 

ht 

stationary above a fixed point on the equator. If this same unperturbed 

satellite is placed in an orbit with an inclination a relative to the 

G 

equatorial plane, it is obvious that the subsatellite point can no longer 
remain fixed since it must vary in latitude between +QC and -a during 
each orbital period. However, if the orbital angular rate is again set 
equal to 0^, the resulting ground trace of the subsatellite point is a 
fixed curve on the rotating earth. This curve is in the form of a figure 
eight with its crossing point on the equator. Thus, although the satel- 
lite itself is no longer stationary relative to the earth, it’s ground 
trace pattern is . 

If perturbations are also considered with the resulting regression 

of the orbital plane, it is necessary to modify the orbital angular rate 

to compensate for this regression. The appropriate orbital angular rate 

can be determined from Fig. 8. This figure shows the intersections with 

the earth's surface of the reference, equatorial and orbital planes. It 

is assumed that the orbital plane is initially in the position indicated 

by the dotted line and that the initial subsatellite point coincides witJa. 

At the end of a time t^ required for n crossings of the reference 

plane, the subsatellite point is at S at an angular distance ijf from 

n ° ^n 

Xj^ . In this same time the point A on the earth, which was the initial 

subsatellite point, is displaced from x by an angular distance 3 . The 

1 n 

quantities t , and 3 can be expressed as follows 
n n n 
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2nTT 

^ e 

o 


(69) 


i|; = ijf t (70) 

n ^ o n 


3 = e„t - 2nn (71) 

n E n 

where 0 is the orbital angular rate to be determined. If the orbital 
o 

angular rate is adjusted so that 6^ is equal to then the subsatellite 

point and A will be coincident whenever the orbital plane passes through 
the axis (i.e., at intervals of half the regression period). The con- 

dition that equals can be expressed by combining Eqs . (69) through 
(71) to give 

e = L - i (72) 

O E O 

With this value of orbital angular rate the resulting ground trace is 
again a figure eight. However, as will be seen, its size and its posi- 
tion relative to the rotating earth are no longer constant. Although 
neither the satellite nor its ground trace are stationary relative to 
the earth, synchronism is maintained since the motion of the ground 
trace pattern in longitude has no secular terms. Thus, Eq. (72) repre— ~ 
sents the condition for synchronism in this case. 

From Eq. (72), it is seen that the desired synchronous orbital 
angular rate depends on the orbital inclination through ijr^. However, 
this variation in 9^ corresponds to a change of the order of ,2 mi 
in orbital radius. This change in r^ is negligible insofar as the 
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determination of either a, or ilr is concerned. 

1 ^o 

Determination of the Ground Trace 

The geometry for the determination of the satellite ground trace 
is shown in Fig. 9, where the angle ^ is the inertial longitude of the 
satellite measured from the axis and Y is the satellite latitude. 

The unit vector i along the satellite radius vector can be ex- 
pressed relative to the x , y , z coordinate system in the form 


i = (cos Y c:os 3) i^ + (cos y sin 3) j + (sin y) k (73) 

G G G 


The same unit vector can be expressed relative to the x^, y^, z^ system 
as 


i = ai-+bj. + ck, 
X 1 X 1 X 1 


(74) 


where a , b and c are given in Appendix A. 

XX X 

If Eq. (74) is projected into the x , y , z system it becomes 

G G G 


1 = a 

X G 


+ l^b^ cos(X - a^) + c^ sin(X - a^)j 


+ 1^- b^ sin(A - a^) + c^ cos(4 - 


( 75 .) 


By equating components of Eqs . (73) and (75), the following ex- 


pressions for 3 ^rid Y are obtained. 


b cos(X - O'.) + c sin(A. - or,) 
X l x 1 


a 

X 


tan 3 ~ 


(76) 
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sin Y = - b sin(>L - a ) + c cos(X - ) (77) 

X XX L 

In Fig. 10 it is assumed that initially the orbital plane passes 
through and the satellite is at x^. The figure then represents the 
geometry at some later time, T, corresponding to some integral number 
of equatorial crossings. S represents the satellite position and A 
is the current position of the initial subsatellite point. The value 
of 9^ can be expressed in terms of by setting y equal to zero and 
replacing \jr and 0 by and 9^ in Eq. (77). Solution of the resulting 
expression gives 


tan 0 

o 


sin \|f sin(A. - 

o 1 


sin O' cos(X - O',) - cos Oi cos \|/ 
o 1 o ’^o 


sin(A. - a^) 


(78) 


where 


ijr = )lf T 


(79) 


From Fig. 10, it is seen that the reference plane crossing im- 
mediately prior to T occurs at time T - t^ , where 


t = 
o 


(80) 


and 9^ is given by Eq. (72). 

If t is the time elapsed since this reference plane crossing, 
then the general expression for i|f in Fig. 9 is given by 


= i (T - t + t) 
^ ^ o o 


( 81 ) 
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The angle 0 in Fig. 9 can also be expressed as 


9=0 (T-t +t)- 2nTT 
o o 


( 82 ) 


where n is the number of orbits completed at the time of reference 
plane crossing. Since n is given by 


n = 


9 (T - t ) 
o o 

2tt 


(83) 


Eq. (82) can be expressed as 


0 = 9 t (84) 

o 

If 0^ represents the current inertial longitude of the initial sub- 
satellite point, A, it can be expressed as 

3 = 0 (T - t + t) - 2nn (85) 

1 JC4 O 

Elimination of 0^ and n between Eqs. (72), (83) and (85) gives 

- t^) (86) 

If Eqs. (72), (81), (84) and (86) are combined, the following value 
for 0^ is obtained : 

0^ = 0 + \|r (87) 

The longitude of the satellite relative to the reference point A on the 
rotating earth is given by 
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A3 = P - 3^ (88) 

while its latitude is equal to y The relative longitude, AP, and the 
latitude, Y, can be computed as a function of t by means of Eqs. (76)- 
(81), (84), (87) and (88). 

For a given value of T, a determination of AP and Y versus t over 
one orbital period determines the ground trace on the rotating earth. 

By taking various values of T over one regression period, the effect 
of the phase of the regression can be shown for various initial orbital 
inclinations as follows. 

Orbit in Reference Plane . In this case, the determination of AP 

and Y be simplified considerably since <y is zero and t is eaual 

o o 

to T. Thus, Eqs. (76), (77), (87) and (88) become 


tan 3 = cos(X - O',) tan 0 t 
1 E 

sin Y = - sin (A - o- ) sin 9 t 



AP = P - Pj^ 


(89) 

(90) 

(91) 

(92) 


The resulting plot of y versus AP is shown in Fig. 11, where it is seen 
that the maximum value of latitude is equal to X - or ±7°20' , which 
is the inclination of the orbital plane to the equatorial plane. Since 
the time T does not appear in Eqs. (89) through (92), the ground trace 
as shown will contine to repeat since this particular orbit remains fixed 
in inertial space. 
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Fig.l 1 — Ground trace (orbit in reference plane, 0) 
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Orbital Plane Normal to Reference Plane . In this case, the equa- 
tions for the determination of AP and Y ^^e simplified by setting 
equal to 90°. Under these conditions, the ground trace is determined 
by the following relations: 

sin(X - O',) sin 0 t + cos(X - O'.) sin \|; cos 0 t 

tan P = i 2 E_ 

COS t cos 0„t 
o E 

sin V = - sin(X - a,) sin ijf cos 0„t + cos(\ - a,) sin 0 t (94) 

i O E IE 

^1 = ♦o 


tan 0 

o 


A3 


sin 


3 - 


tan(X - Ck'^) 


(96) 


(97) 


As shown in Appendix F, may oscillate about a value of either 90° 
or 270°, with an amplitude less than 90° and a period greater than 
269 years . 

By means of these equations, the ground traces shown in Fig. 12 

have been computed for equal to 0°, ±90° and 180°. It is seen that 

o 

II o o 

as ItQl increases from 0 to 180 the amplitude of the figure eight in- 
creases from 90 - (A. - a^) , or 82 40' , to 90 ± (X - or^) , or 102 40' . At 
the same time the maximum latitude ranges from 82°40' at equal to 0° , 
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to 90° when equals 90°, and back to 82°40' for equal to 180°. 

The longitude variations for < 90° are oscillatory in nature, but, 

for > 90°, the longitude A3 increases monotonically . This is due 

to the fact that the ground trace encircles the earth's axis in a 
negative sense for these latter conditions. 

It should be noted that only the ground trace for equal to 90° 
is actually fixed relative to the earth. Those corresponding to other 
initial values of are subject to large amplitude oscillations in 
longitude . 

Orbit in Equatorial Plane . In this case the orbital plane is 

initially in the earth's equatorial plane, with set equal to A. - 

which is equal to 7°20' . From Eq. (59), it can be shown that the or- 

-4 

bital regression rate is equal to -3.231 X 10 rad/solar day, which 
corresponds to a regression period of 53.249 years. The ground trace 
can now be determined by means of Eqs . (76) through (78). Since the 
trace changes as the orbit regresses, the computation is made at five- 
year intervals in T up to 25 years. The resulting patterns are presented 
in Fig. 13, where it is seen that during the first half of the regression 
period the ground trace grows from a single point to a figure-eight pat- 
tern which attains its maximum size after half of one regression period. 

At this time, the ground trace has dimensions of ±14°40' in latitude and 
±1° in longitude. In addition to the variation in size of the earth trace, 
its equatorial crossing moves relative to the origin, which is the initial 
subsatellite point. During the first half of the regression period, the 
equatorial crossing moves to the east, reaching a maximum displacement 
of about .6°, which then decreases to zero after half the regression 
period. The behavior of the ground trace during the second half of the 
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Fig. 13 — Ground trace (orbital plane initially equatorial. 
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regression period is the reverse of that just described in that the 

pattern decreases in size and degenerates to a point at the end of one 

regression period. Also, the equatorial crossing drifts toward the west, 

o o 

reaching a maximum displacement of .6 and returning to 0 at the end 
of the regression period. The behavior of the ground trace during the 
second half of the regression cycle is a mirror image of that shown in 
Fig. 13. 

If the geometrical representation used in Fig. 6 is adapted to this 

case, the relationship between the z, z and z axes and the trace of the 

J. G 

z axis on the reference sphere is shown in Fig. 14. Since and A. - 

are equal, the earth's polar axis intersects the z axis trace. Also, the 

arc of the great circle through z and z is a measure of the instantaneous 

G 

inclination. Of , of the orbit relative to the earth's equatorial plane. 

G 

Thus, as the z axis moves around its trace, the value of varies from 

0^ at z to a maximum of 2(A. - ) or 14°40’ at a point diametrically 

G JL 

opposite to z . This instantaneous value of O' is equal to the maximum 
G G 

latitude excursion of the ground traces of Fig. 13. 

Orbit at an Arbitrary Inclination . As an example of the general 

behavior of an inclined synchronous orbit, a value of 30° for is 

selected. This inclination corresponds to a regression rate, ij/^, of 
-4 

-2.821 X 10 rad/solar day, which results in a regression period of ^ 
60.98 years. The resulting earth traces at five-year intervals up to 30 
years are shown in Fig. 15. It is seen that initially the characteristic 
figure-eight pattern has an amplitude in latitude of - (A - or^) or 
22°40'. During the first half of the regression period, this amplitude 
increases to a maximum equal to + (X - or 37°20' . During this time 
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Fig. 14 — Regression of a synchronous equatorial satellite 
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Fig ,15 — Ground trace ( a 
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the equatorial crossing point, which is initially at the origin, drifts 
in an easterly direction, reaching a maximum of about 2 , after which the 
drift reverses and the crossing point returns to 0 at the end of one- 
half of the regression period. As in the previous case, the behavior of 
the ground trace during the second half of the regression period is a 
mirror image of that during the first half. 

The relationship of the z, and z^ axes is shown in Fig. 16 for 
this case. It is seen that the instantaneous orbital inclination, 
relative to the equatorial plane lies between - (X - and 

O' + (X - Q'^) , which are also the limits of the maximum daily latitude 
excursions of the ground traces of Fig. 15. In general, this relation- 
ship can be stated as 

|q'^ - (X - Ck'^) I ^ ^ + (X - cy^) (98) 

Size and Position of the Ground Trace 

It is seen from the previous examples that the size and position 
of the ground trace are functions of both the orbital inclination, > 
and the time. In order to describe these variations, it is convenient to 
determine both the maximum latitude of the ground trace and the position 
of the equatorial crossing as functions of and time. 

Variation of Maximum Latitude . Since the maximum latitude during 
a single traversal of the ground trace is equal to the inclination, cy^, 
of the orbit to the equatorial plane, this maximum latitude, 
expressed as 

cos Y = (k • k ) 
m G 


( 99 ) 
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which can be expressed from Appendix A as 


cos y == cos O' cos(X - o',) 
% o 1 


+ sin a sin(X - cv, ) cos i|/ 
o 1 ' 


( 100 ) 


where 


= ilf T 
^o 


( 101 ) 


and is given by Eq. (59) . 

The dependence of y on T and a is shown in Fig. 17, where T has 

m o <=> j 

been normalized by dividing by T . It is seen that the maximum possible 

K 

variation of y for a given value of O' is 14^40' . In addition, for 
m o 

those contours for O'^ > 82 40’, the central portion corresponds to orbits 

which encircle the earth's pole in a negative direction. 

A further examination of Fig. 17 shows that, with the exception of 

the contour for (x = A, - ck' , all of the o' contours have a zero slope at 
o 1 o 

T — 0. However, for = X - o(^ the slope is .8630*^ per year, which also 
represents the initial rate of change of the orbital inclination relative 
to the equatorial plane. 

Figure 18 shows both this variation of y for O' equal to A - o', , 

m o 1 

and the equivalent curve from Fig. 10 in Ref. 1. It is seen that the ^ 
present analysis gives a value of 53.249 years for the period of the lati- 
tude variation and a maximum amplitude of 14°40' , compared with the values 
of 73.6 years and 20° obtained by the more approximate methods of Ref. 1. 
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Fig. 17 — Maximum latitude as a function of and T 
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Equatorial Crossing Position . The variation in the position of the 
equatorial crossing point as a function of time and orbital inclination 
can be determined as follows. At the time of equatorial crossing, the 
inertial longitude of the satellite can be obtained from Eq. (76) in the 
form 


cos 9 sin \lj cos(X-Q', ) 
o ^ 0^1 


+ s in 


tan 3 = 

o 


in O' sin(X-Q'.) + cos a cos (A.-O', )cos il; 

q 1 q ^ 1 ^o_LI 


cos 0 cos lif - sin 9 cos o' sin ik 
o o o o o 


( 102 ) 


Elimination of 0^ between Eqs. (78) and (102) gives 


tan 3 


sin ilf sin O' 
o o 


o - cos a sin(A. - o’.) + sin ot cos(X - cy, )cos \1; 

o ' 1 o ' 1' ^o 


(103) 


The longitude of the initial subsatellite point is given by Eq. (87) 
as 

3, = e + ^ (104) 

lo o ^o 

Thus, the longitude difference between the current equatorial crossing 
and the initial subsatellite point is given by 

A3 = 3 - 0 - (1^) 

o o o ^o 

By means of Eqs. (79), (103) and (105), A3^ can be computed as a function 
of and the time T. The resulting relation is shown in Fig. 19, where 
T has been normalized as in Fig. 17. It is seen that the resulting sur- 
face is such that for a given value of the quantity has a maximum 

at approximately one-quarter of the regression period and a minimum at 
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Fig. 19 — Longitude of equatorial crossing as a function of and T 
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about three fourths of this period. The value of this maximum ranges 
from zero at equal to 0° to approximately X - at equal to 
90°. The exact expression for this maximum and its time of occurrence 
can be found by differentiating Eq. (105) with respect to time and 
equating to zero. This results in 




= tan 


-1 


max 


Of 

Q ■ O . 

o sin -j- sin 


iV-Qr s /cos a + cos (A-a ) 


[3 cos(X-d ) + 3 cos O' - cos a cosCA-o-, ) - l] 
1 o o 1 


(106) 


and 


max _ 1 -If “o 

T~ ~ 27F L’ ~ 

R 


A - O' 


^)J 


(107) 


It is seen from Eq. (107) that the maximum occurs very slightly past 
one-quarter of the regression period. In a similar manner it can be 
shown that the minimum of has a magnitude equal to the negative of 

Eq. (106) and a time of occurence given by Eq. (107) as slightly before 
three-quarters of the regression period. Contours for these maxima and 
minima are shown in Fig. 19. 

It should be noted that Fig. 19 applies only to regression about 
the axis and that the effect would actually be superposed on the 
shorter period oscillation due to equatorial ellipticity described in 
Ref. 1. 

Period Between Equatorial Crossings 

The nodical period is defined as the time interval between crossings 
of the reference plane at the ascending node and is given by 
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T 

O 


2tt 2it 


0 

o 



(108) 


However, the period between successive equatorial crossings is not 

necessarily equal to when orbital regression is present. The amount 

of this deviation from t can be determined as follows. Figure 20 shows 

o 

the position of the orbital plane at the times of two successive crossings 

of the equatorial plane. If these two crossing times are T and 

T + T + At, then the associated values of \lf are \|f and ilr + Ailf , and 
o o o o 

those for 0 are 0 and 0 4- A0 . The value of the increment Ai is given 

o o o o 

by 


Aijr = \jr (t 4- At) 
o o o 


(109) 


On the other hand, since the angle 9 increases by approximately 2n be- 
tween two equatorial crossings, the angle A0^ can be determined from the 
relation 

0 4- a9 =0 4- 0 (T 4- At) - 2n (110) 

oooo o 

which together with Eq. (108) gives 


A0 = 0 At (111) 

o o ^ 

At the time of the first equatorial crossing, T, the latitude is 0°. 

Thus, substition of ilr and 0 for il; and 0 in Eq. (77) gives 

o o 


sin(X - O',) cos 0 sin \1; 4- sin 0 cos 

^ 1 L o o o 


4- sin 0 sin Ck* cos(X - ) = 0 

o o 1 


llf cos 
^o 



( 112 ) 
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Similarly, at the time, T + + Z\t, of the next equatorial crossing, 

the latitude is again 0°, and the substitution of ilr + Aijf and 0 + A0 

o o o o 

for and 0 in Eq. (77) gives a second relation of the form 

- sin(X - O', ) cos(0 + A0 ) sin(\l; + Ailf ) 

1 L o o ^o ^o 

+ sin(0 + A0 ) cos(t|r + Aijf ) cos 1 

o o o o oj 

+ sin(0^ + ^o ^ (113) 


By expanding Eq. (113) for small values of A0^ and Aijf^ and subtracting 

Eq, (112), the following relation between A0 and A^ is obtained, 

o ^o 


A0 (sin 0 sin ijf - cos 0 cos ijf cos cx ) sin(A - ) 

o L o o o o 0 ^1^ 

+ cos 0 sin O' cos(A - O', ) 
o o' I'J 



cos 


0 cos llf 
o ^o 


sin 


0 sin il; cos o' sin(A - oc. ) = 0 
o ’^o oJ 1 


(114) 


Elimination of A0^ and A^^ between Eqs . (109), (111) and (114) gives the 
following solution for At. 


At = 



T sin(X-cy,) 
o '1 

- cos O' sin(X-Q',) + sin a cos(X- 0 ',) cos \|f 
L o ' 1 o ' 1 ^oj 


0 

o 

r 2 2 

1 sin sin 

(X-Ck'j^) + sin cx^ cos(X-a^) | 


- cos cx^ sin(X-Qrj^) cos 


I + sin(X-a^) cos <x^ sin(\-Qi'^) - sin cos(X-cy^) cos 


( 115 ) 
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where \lf^ is given by Eq. (79). 

Figure 21 shows the relationship between At, and T specified by 
Eq. (115). It is seen that the period between equatorial crossings dif- 
fers from T by only a few seconds for most values of Ck^ and T. However, 
o o 

in the vicinity of the singularity at = A. - and T = 0, At can be 
appreciable since under these conditions the orbital plane and the equa- 
torial plane are very nearly coplanar and the position of equatorial 
crossing can shift very rapidly. 

ORBITAL INCLINATION CONTROL 

In the preceding discussion it has been shown that in the absence of 
any control forces the characteristic figure-eight ground trace may have 
long period variations in size and position due to orbital regression. 

For many applications it may be desirable to limit these variations. 

This is particularly true of the excursion in latitude where, for coverage 
reasons, it may be necessary that the satellite latitude excursion be re- 
stricted. This sort of control can be accomplished in two ways, the 
most obvious being to apply control forces in such a way as to reduce 
the orbital regression rate to zero, with the result that the orbital 
plane remains fixed in inertial space and the ground trace remains fixed 
in size and position on the earth. The other method of control is passive 
and is achieved by injecting the satellite into orbit in such a way that 
even in the presence of regression the period during which the latitude 
remains less than a given value is maximized. 
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Fig. 21 — Period between equatorial crossings as a function of and T 
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Active Control 

In connection with this method it is of interest to determine the 
amount of propulsion required to reduce the regression rate to zero 
thus fixing the ground trace relative to the earth. 

If it is assumed that a control force F is applied in the z 
direction perpendicular to the orbital plane, Eq. (21) has an additive 
term of the form tight side, and the resulting expression 

for the regression rate corresponding to Eq. (48) becomes 


SJ^R^e 


t = - 


+ 


2 

sin O’ 


r 


o 

o 






30 



m 

(■ 



M-0 

sin O' 


o o 



a2 



30 

(r 

• 


e 

sin cy 


o 

o 


F 

sin 0 


cz 



m 


i) (i * k) sin 0 
m 


M 9 sin a 
so o 


(116) 


The resulting steady-state regression in the presence of the force^ 

F is given by 
cz 

T 

♦ c = I J ♦ dt 017) 

o 

where T is a time large compared to any of the oscillatory periods 
associated with ijr. Combination of Eqs . (116) and (117) gives 
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ij; + 

^ o 


r M 9 T sin O' o 
o s o o 


F sin 9 dt 
cz 


(118) 


To reduce the steady-state regression rate, to zero, it is necessary 

that 




a sin 9 dt 
cz 


- r 9 Tlf sin O' 
0 0^0 o 


(119) 


where a^^ is the control acceleration perpendicular to the orbital plane 

From the form of Eq. (119) , it is seen that accelerations applied 

when 9 is 90° are more effective in controlling the regression rate. 

Thus, it is assumed that the control consists of a single impulsive 

o 

force applied each time 9 is equal to 90 . Under these conditions the 

left side of Eq. (119) represents the required velocity change per unit 

time, AV /T, so that 
z 


AV 


-rr- = - r 9 it sin a 
T o 0^0 o 


( 120 ) 


Combination of Eqs. (59) and (120) gives 


AV 3r 0 sin 2a 
z o o 

T 16 


fl H (2 -3 sin o')l( 2 - 3 sin or.) 

L o m J i 

2[ji0 


‘2 2 
2J^e R 




J V lx r" o “ 

, 22 " ^ ■ ^ 1 ^ 


0 r 


(121) 
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Thus, the propulsion requirement as specified by Eq. (121) depends only 
on the inclination angle, relative to the reference plane. 

Since the orbital inclination angle, Oi . relative to the equatorial 

Cj 

plane determines the maximum latitude excursion of the satellite, is 

the parameter which would be specified. The geometrical relationship 

between and or is shown in Fig. 22, where it is seen that or varies 

depending on the position of the z axis on the constant ot contour. 

G 

However, 0(^ ranges between extremal values of {X - and 

(X - - ag. These values occur when ^ is equal to 0° and 180° re- 

spectively. Thus, the propulsion requirement is also bounded by the 
expressions obtained by substituting the above extremal values for 
in Eq. (121). The propulsion requirement can then be expressed as 


AV 

sin 2(X - ± Og) (122) 


where is the amplitude of Eq. (121) and the use of the plus or minus 
sign is determined by which of the two gives the smaller absolute value 
for AV /T. The propulsion requirement as a function of is shown in 
Fig. 23, where the solid lines represent the minimum condition speci- 
fied . 


It is seen that for a equal to zero (equatorial orbit) the propul»- 

G 

sion requirement is 151.9 ft/sec/year, which is in good agreement with 
the value of 150 ft/sec/year obtained in Ref. 1. For a equal to 7°20', 
the orbit lies in the reference plane and remains there by definition. 
Thus, there is no propulsion required for orbital control. As cy_ in- 
creases, the propulsion requirement reaches a maximum of 580.4 ft/sec/year 


when 


equals 45 . Although Fig. 23 indicates a zero value for the 
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Fig. 23 — Propulsion requirement for orbit inclination control 
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propulsion requirement when a equals 82°40' , this orientation is of 
little interest since the orbit is in unstable equilibrium, as shown in 
Appendix F. However, on the basis of Appendix F, there is a second 
stable orbital orientation with the orbital plane perpendicular to the 
direction of the vernal equinox. This orbit also has a zero propulsion 
requirement . 


Passive Control 

In this method of ground trace control it is necessary to select an 

Initial phase of the regression cycle such that the latitude excursion 

is bounded for a maximum period of time. To accomplish this, use is 

again made of the reference sphere as shown in Fig. 24, where the circle 

about z_ is a locus of constant a and thus constant y as defined in 

G m 

Eq. (100). The circle about is the path of- the z axis as the orbit 
regresses. If the z axis is initially at A, then the maximum latitude 
excursion will remain less than its initial value, V » during the time 
it takes the z axis to regress from A to B. This time can be obtained 
by solving Eq. (100) for T to give 


T 



cos V " cos <y cos(X - O',) 
^mo p_ 1 


sin a sin(\ - o',) 
o 1 


(123) 


By differentiation of Eq. (123) it can be shown that for a given value 

of Y the maximum value of T occurs when a is given by 
'mo o ® 


cos cy = 
o 


cos(\ - cv^) 


( 124 ) 




Fig. 24 — Passive orbital control 
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In obtaining Eq. (124) the dependence of on (y^ is neglected since 
only values of q>^ less than are considered. The resulting ex- 

pression for the maximum value of T becomes 


\l2 sin(\ - a. + Y ) sin(X - o'. ) 

1 I—IeS ^ 1 S2_ (125) 

sin(A - 

In Fig. 25, is plotted as a function of the maximum permissible lati- 
tude, V • It is seen that for or equal to A. - o', , T reaches a value 
^mo o ^ 1 m 

of 26.4 years . 

For values of v greater than X - o', , T is infinite since the 
'mo i m 

7 contour of Fig. 24 encloses the z axis, which is one of the naturally 
stable positions for the z axis. 

Thus, for any given satellite application, if the specifications 
of maximum permissible latitude variation and minimum acceptable life- 
time are represented by a point to the left of the curve in Fig. 25, 
then it is necessary to use an active control system with its associ- 
ated propulsion requirements as shown in Fig. 23. 
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IV . CONCLUSIONS 


As a result of the analysis presented in this Report, the follow- 
ing general conclusions can be stated regarding the effects of the sun, 
the moon and the oblateness of the earth in perturbing satellites in 
nominally circular orbits. 

o For a given orbital radius, an orientation of the orbital 

plane can be found which remains essentially constant relative 
to inertial space under the influence of the assumed per- 
turbing forces . 

o The plane of the ecliptic, the earth's equatorial plane and 
the constant reference plane defined above have a common line 
of nodes, and the inclination of the reference plane relative 
to the ecliptic is less than that of the equatorial plane, 
the exact value being a function of orbital radius. 

o If an orbit is inclined relative to its reference plane at 
an angle less than 79°, its orbital plane will move in such 
a manner that its inclination relative to the reference plane 
remains essentially constant while its line of nodes in the 
reference plane regresses at a steady-state angular rate pro- 
portional to the cosine of the inclination angle. 

o A highly elliptical type of regression can take place about 
the direction of the vernal equinox if the inclination of the 
orbit exceeds 79°. 

o In addition to the steady-state motion, there are a large num- 
ber of oscillatory terms in both inclination angle and re- 
gression angle. These result in an oscillation in the direction 
of the normal to the orbital plane of less than .5° relative 
to its steady-state motion. 

In addition to these general conclusions, which apply to orbits^of 
any radius, certain others are reached which apply only to synchronous 
altitude orbits regressing about the normal to the reference plane 

o The reference plane for a synchronous altitude orbit has an 
inclination of 16°7' relative to the ecliptic and -7°20' rela- 
tive to the equatorial plane. 

o The regression period of a synchronous orbit is 52.81 years at 
an inclination of 0^ and varies inversely as the cosine of the 
inclination angle, relative to the reference plane. 
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o The ground trace of an orbit lying in the reference plane is 

a figure eight with crossing point on the equator and a maxi- 

mum latitude excursion of 7°20' . This trace is invariant in 
size and location as a function of time. 

o An orbit inclined to the reference plane has a ground trace 

which is also a figure eight, but both its maximum latitude 
excursion and its equatorial crossing point undergo oscillatory 
variations with a period equal to the regression period. 

o In the case of a synchronous orbit which is initially equatorial, 
it is found that the ground trace develops from a point to a 
figure eight with a latitude excursion of 14®40' after 26.6 years 
and then shrinks back to a point after another 26.6 years. 

o To stop the orbital regression and thus fix the satellite ground 
trace requires a fuel expenditure proportional to the sine of 
twice the inclination angle relative to the reference plane. 
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Appendix A 
DIRECTION COSINES 



The following tables present the direction cosines relating the 
unit vectors of the various coordinate systems defined in the body of 
the report . 


REFERENCE COORDINATE SYSTEM 



i 

o 

L 

k 

o 


1 

0 

0 




sin Q'^ 


0 

cos 






0 

- sin 

cos on^ 


EQUATORIAL COORDINATE SYSTEM 





k 

1 


1 

0 

0 


0 

cos(X - Q'^) 

sin(\ - 


— 

0 

- sin(X - 

cos(X ' 
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SATELLITE ORBITAL SYSTEM 



LUNAR ORBITAL SYSTEM 
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Appendix B 
EQUATIONS OF MOTION 


In the body of this Report, the equations of motion of 
the earth and the moon relative to inertial space are given 
(15) and (16) in the form 


the satellite, 
in Eqs . (14) , 


Sate llite : 


..FEE 
“ _ S s ms Es 

^s M M M 
s s s 


(B-l) 


Earth: 




(B-2) 


Moon : 


- _ ^Sm , ^Em 

R = z: ^ z : — 

m M M 
m m 


(B-3) 


The vector equation of motion of the satellite relative to the earth is 
obtained as the difference between Eqs, (B-l) and (B-2) as given in 
Eq. (17) and below 


r = 


F_ + F + F^ 

S s ms ^ 

M 


F + F 
SE mE 




(B-4) 


Similarly, the vector equation of motion of the moon relative to the ear-th 
is given by the difference between Eqs. (B-2) and (B-3) as 

F + F F ~ 

- ^ Sm Em SE 

^o M " 

m 

Finally, the vector equation of motion of the center of mass of the earth- 
moon system relative to inertial space is obtained as follows. The 


+ F 


mE 




(B-5) 
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position of the center of mass is located by the vector R, which is de- 
fined by the relation 


(Me + M^) r = MeRe + M^R^ (B-6) 

Elimination of and R^ between Eqs. (B-2) , (B-3) and (B-6) gives the 
desired equation for center of mass motion as 


R 


F + F + F + F 
SE Sm Em mE 




+ M 


(B-7) 


GRAVITATIONAL ATTRACTIONS 

The various gravitational forces involved in the foregoing equa- 
tions are evaluated as follows. 


Force of the Sun on the Satellite 


The force F„ can be expressed as 
S s 


GM^M 

~ b s ~ 

F = ' T 

Ss 3 S 


where 


r = R - — + r 

S M- 


The magnitude of r is then given to first order as 


rg = + 2(r . R) . 2(-^ • R 


(B-8) 


(B-^) 


(B-10) 
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from which the inverse cube term can be expressed by a binomial ex- 
pansion as 



3 — — 3 

- ^ (r • R) + — 

R R 




(B-11) 


Substitution of Eqs . (B-9) and (B-11) in Eq. (B-8) gives 



GM^M 

h S 


R 





+ r 



R) R 



r) r] 


(B-12) 


neglecting terms of the order of r /R . 


Force of the Moon on the Satellite 


The force F is given by 
ms 


GM M 
m s 

3 

r 

m 


where 


F = - 
ms 


(B-13) 


r = - p + r 
m o 


(b-147 


From these relations the first-order expansion 
the same manner as above in the form 


for F 

ms 


is obtained in 


F 

ms 


GM M 
m s 


p + r 
o 


3(P. 



3 


2 


(B-14) 
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2 2 

neglecting terms of the order of r / • 


Force of the Earth on the Satellite 

The gravitational potential of the earth at a given point can be 

expressed in terms of the radial distance r from the earth’s center and 

the distance z from the equatorial plane as follows 
G 


u 







0 


(B-15) 


where G is the universal gravitational constant, is the coefficient 

-3 

due to earth oblateness (J 2 = 1.08219 x 10 ) and is the mean radius 

of the earth. The force components in the r and z directions are found 
by differentiation to be 


and 


F 

r 







2r‘ 




F 

z 



5 

r 


S ince 


z 


G 


(r • 




(B-16) 


(B-17) 


(B-18) 
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the force exerted by the earth on the satellite can be expressed in 
vector form by combining Eqs. (B-16) through (B-18) to give 


GMgM^ 


Es 


5(r . k ) 

1 - ?T- - 1 r 




3J„R 

. 1 Q — 

+ — ^ — (r 


k ) k 
G G 


(B-19) 


Force of the Sun on the Earth 


The force F can be expressed as 
bE 


SE 




(B-20) 


where 




= R - 


(B-21) 


Again, by a binomial expansion, F becomes 

SE 




SE 


R' 


__ P ^ _ _ _ 

R + 7T (p • R) R 


(B-22) 


2 2 

neglecting terms of the order of p^/|J.R . 


Force of the Moon on the Earth 


In this case the force F^ is given by 


m£ 


^m^ - 
3 Po 


(B-23) 
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Force of the Sun on the Moon 


The force F„ is given by the relation 
Sm 


Sm 


R 

J m 


R 


m 


(B-24) 


where 


R = R + (l - p 
m \ U/ o 


(B-25) 


These two relations can be expanded as before to give 


Sm 


GM^M 
S m 


— / IN — 

+ 1' - po 


R 


3 

T 


- B (po • 


(B-26) 


2 2 

neglecting terms of the order of P^/R • 


Force of the Earth on the Moon 


In this case the force F^ is given by the relation 

Em 


GM^M 

Em-* 


Em 3 ^o 


(B-27) 


MOTION OF THE MOON 

Substitution for the forces on the right side of Eq. (B-5) gives 
the following equation for the motion of the moon relative to the earth: 
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G(M^ + - 

3 Po 

Po 

GM r_ •? _ _ -n 

--^L^o -i(Po • 

R R 


which can be simplified to 


Po = - 


G(Me + M^) _ 


(B-29) 


since the second term in Eq. (B-28) is less than one percent of the first 
term. 


Since it has been assumed that the moon rotates at a fixed angular 
rate and at a constant distance from the earth, its acceleration is given 
by 


11 '*2 — 
p = - 9 p 

o mo 


(B-30) 


which together with Eq. (B-29) gives 


where 



m 









+ M 


m 


M 

m 


(B-31) 


(B-32) 


MOTION OF THE EARTH -MOON SYSTEM 

In a similar manner, the forces in Eq. (B-7) can be eliminated to 


give 
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R 


(B-33) 


as the equation of motion for the center of mass of the earth-moon system. 
Since this motion is also assumed to have a uniform angular rate, and 
a fixed value of R, the acceleration is given by 



(B-34) 


which combined with Eq. (B-33) gives the following for the earth's or- 
bital angular rate 


©' 



(B-35) 


MOTION OF THE SATELLITE 


Vector Equation of Motion 

Substitution of the appropriate force expressions in Eq. (B-4) and 

elimination of the quantities GM and GM by mean of Eqs. (B-31) and 

o in 

(B-35) give the following vector equation of motion for the satellite: 



(B-36) 
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which corresponds to Eq. (18). 

Substitution of ri, Rr. and p i for r. R and p respectively in 

1 o m '^o 

Eq. (B-36) gives 


r = - 


GM^ 

2 




3J 

2 o 

n 2 

2r 


5(1 • 


-.j)t 


3J,R 

+ — ‘'g 


- r0 


9^[l - 3(1 • r^) 


r0' 


m 


y ■ • V 


(B-37) 


The left side of Eq. (B-37) can also be expressed in the form 


j2— r- 


uo X X r] 
o o J 


(B-38) 


where the derivatives on the right are relative to the x, y, z coordinate 
system, and (ju^ is the angular velocity of the x, y, z system relative to 
inertial space. 

Substitution of ri for r in Eq. (B-38) gives 


7 = ^ I + 2 ^ X T 
dt2 dt L O . 


+ r 


.d(D 

o T 1 
-T- X 1 

dt J 


+ rl ((ju • i) cju - (uu • (ju ) - , 
L o o o o J 


Tl 


(B-39) 
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Component Equation of Motion 

If the i, j and k components of Eqs. (37) and (39) are equated, 
the following equations of motion are obtained 


— 7T + r (uo 
dt^ L o 


T)^ 




3GM^J2R; 


2r 



- r^' 



r0 


H r 1 



(B-40) 


2 


dt 



k) + r 




+ r(o) 
o 


i) (o)^ . j) = 


4 



i) • j) 


+ 3r@^(r^ • i) (r^ • j) 


3r0' 


m 


(i 


m 


i) (i 


m 


j) 


+ 


M- 


(B-41) 
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H — — — \ _ _ _ _ 

2 ^ - KdT ■ j) + ■ '^'> 


,r 


4 <*^0 ■ 


+ 3r@^(r^ • i) (r^ • k) 


3r9^ _ _ _ _ 


(B-42) 


These three equations can be further simplified to 


dr r.“- 

dt^ ' 


— — 2 — “ 2”! 

• j) + (% • ^) J = - — 


3GM^J2 


R 


2r 


4 L 


T - 3(k^ • 1)2; 




r 

I 1 
L 



(B-43) 
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1 d r 2 

r dt 




+ r(cjD 
o 


i) 


• j) = 


3GM^J2R„ 


1) (kg 


j) 


+ 3r9^(r^ • i) (r^ • j) 


3r9 


m 


(i 


m 


i) (i 


m 


j) 


(B-44) 


_ i -i. r 2 

r dt 





+ r(cjD 
o 


1) («)^ • k) = 


3GM^J2R^ 

4 

r 



i) • k) 


+ 3rf0 (r- • i) (r. 


k) 


+ 



M- 


(i 

m 


i) (i • k) 
m 


(B-45) 


Constraint Equation 

In order that the xy plane always be the instantaneous orbit plane, 
it is necessary that the vectors r and r both lie in this plane. The 
vector r satisfies this condition since it is along the x axis by 
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definition. In the case of r the following relation must be satisfied 

(r . k) = 0 (B-46) 

which states that r is perpendicular to the z axis and thus lies in the 
xy plane as desired. 

Equation (B-46) can be simplified by combining it with the expression 

r = ~ + [o)^ X r] (B-47) 

which gives 

(^ * k) + X r] • k) = 0 (B-48) 

Substitution of ri for r reduces Eq. (B-48) to 

^ (i ♦ k) + r(uj^ • [i X k]) = 0 
or 

(^ • T) = 0 

Complete Equations of Motion 

By means of Eq. (B-50) , the three equations, Eqs, (B-43), (B-44) 
and (B-45), can be simplified to 


(B-49) 


(B-50) 
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d^r i-,2 

^ - r(.^ . U) - - ^ 


3GM^J2Ro 


2r 


[1 - 3(k^ • i)^] 


r@^[l - 3(r^ • i)^] 


re 


r - 3(v • 


(B- 51 ) 


3GM J R ^ _ 

[r"(cr^ • k)] (k^ • i) (k^ • j) 


^ - 2 .- 
dt " ' o 


+ 3r^e)^ (r^ • i) (r^, • j) 


2 • 2 

3r 0 __ 

H — (i • i) (i • j) 

|jL m m 


(B-52) 


(o)^ . i) (o)^ . k) = - (k^ • i) (k^ ‘ k) 


+ 30 (r^ • i) (r^ • k) 


30 _ 

+ ^ (i . i) (i . k) 

[J. m m 


(B-53) 
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The three angular velocity components can be expressed by means 
of Eq. (13) so that 


(uu^ ♦ i) = \jf sin 0 sin o’ + cy cos 0 


(cr^ • j) = cos 0sinQ/-Q'Sin0 


• k) = 0 + \]/ cos (y 


(B-54) 

(B-55) 

(B-56) 


Substitution of Eqs . (B-54) through (B-56) into Eqs . (B-50) through 
(B-53) gives 


- r(0 + \jj cos ca)^ 
dt 



2 


r 


sgm^JzR; 


[1 - 3(k, 


1)^] 


2r 


- r@^[l - 3(r^ • i)^] 

_ _ 2 

- [1 - 3(i • i) ] (B-57) 

(i m 


, 2 . 3GM J 

^ [r (0 + ijf cos cy)] = - ^ (k^ • i) (k^ • j) 

r 

+ 3r^0^(7^ ■ I) (7^ . j) 


2*2 
3r 0 

+ (i 

|j. m 


i) (i 


m 


j) 


(B-58) 
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2 .- 


3GMgJ2Ro(kc * i) * R) sin 0 

_ ^ ^ 

r (0 + ijf cos O') sin 


+ 


3©^(r^ • i) (r^ • k) sin 9 
(0 + \jr cos cy) sin a 


+ 


30^(i • i) (i • k) sin 0 

m m m 

M-C© + cos cy) sin cy 


(B-59) 


3GMgJ2R^(k^ • i) (k^ . k) cos 9 
r^(9 + \{r cos cy) 

30^(r^ • i) (r^ • k) cos 0 
(0 + ijr cos (y) 

30^(i • i) (i • k) cos 9 
m m m 

|J. (0 + cos cy) 


(B-60) 


These four equations represent the desired equations of motion, and 
the scalar products of the unit vectors can be evaluated by means of the 
direction cosines given in Appendix A. 
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Appendix C 

OSCILLATORY AMPLITUDES 


INTRODUCTION 

In the main part of this Report, expressions are derived for or, 
the rate of change of orbital inclination, and ij/ , the orbital regression 
rate, in the form 

aCM^J^R^Ck^ • i) (k^ • k) sin 0 
r (9 + ijf cos O') sin O' 


3@^(r^ • i) (r^ . k) sin 0 
(0 + i|f cos O') sin cv 

3^(i • i) (i . k) sin 0 

m m 

|x(0 + ijf cos a) sin (y 


(C-1) 


3GMgJ2RQ(kg • i) (k^ • k) cos 0 

5 * * 

r (0 + \]/ cos O') 

3@^(r^ • i) (r^ • k) cos 0 
(0 + l|; cos Ck') 

30^(i • i) (i_ • k) cos 0 

m m ^ 

|i(0 + \i( cos O') 


(C-2) 


corresponding to Eqs . (37) and (38). It is also indicated in Eqs. (49) 
and (50) that these expressions can be represented as 
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127 

+ = I 

° 1=1 

and 

127 

O' = ^ (C-4) 

i=l 

In this appendix the method of determining the quantities ijj^, , 
and uj^ is indicated and the resulting expressions are tabulated. 

METHOD OF COMPUTATION 

If the scalar products on the right sides of Eqs. (C-1) and (C-2) are 

evaluated by means of the direction cosines in Appendix A, the resulting 

expressions can be expanded as a constant plus a summation of cosine terms 

in the case of Eq. (C-1) and as a summation of sine terms in Eq. (C-2). 

The arguments of these sine and cosine terms are in the form of 

linear combinations of the angles 0, 0, 9 , i and ilr , while their ampli- 

m m 

tudes are functions of the angles o', a and o', as well as the coefficients 

m 1 

of the terms in Eqs. (C-1) and (C-2). If it is assumed that the oscillatory 
variations of ij; are small compared to its steady-state value, then 

the angle \j; on the right sides of Eqs. (C-1) and (C-2) can be replaced- 
by Similarly, if it is assumed that the oscillatory variations of 

0 ( are small, then or can be replaced by its steady-state value, 0(^, 

Finally, it is assumed that t is small compared to 0 so that the 
term 9 + \jr cos Qf can be replaced by 0^. Since the solution of Eqs. (C-1) 
and (C-2) for ijf and Oi depends upon these assumptions, it is essential 
that the resulting solutions verify the assumptions; if this were not 
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the case, the method would be invalid. This aspect of the problem is 
investigated in Appendix D. 

If, in addition, the following substitutions are made 


0 = 


0 t 
o 


0 = ©t 


0 = 0 t 

m m 


'If 


m 


tI; t 
^m 


(C-5) 


then Eqs. (C-1) and (C-2) take on the form indicated in Eqs. (C-3) and 

(C-4) , where the uu. values are linear combinations of the angular rates 0 , 

1 o 

0 , 0 , )tr and i|/ , and the quantities \k , A. and B, are functions of a , 
m ^m ^o ^ ^o 11 o 

(y and a. . 
m 1 

The right sides of Eqs. (C-1) and (C-2) each consist of three terms 
arising from the effects of earth oblateness, the sun and the moon. 

Thus, in the determination of \k as well as A. and B, , the presence of 

^O 11 

J- indicates an oblateness effect, while 0 and 0 indicate solar and 
/ m 

lunar effects respectively. Actually, as will be seen, and some 

of the A.’s and B. 's have contributions from all three. 

1 1 

ANALYTICAL EXPRESSIONS 

Since the details of the actual analysis are rather lengthy, only 
the resulting amplitudes and frequencies are given here. In presenting 
these expressions it is convenient to make the following definitions: 



96 


K = 


o 


(C-6) 


P = 


3® 

89 


1 + (2 - 
2(J,0 


Q • 2 . 

3 s in O' ) 
m j 


(2-3 sin Q'^) 


+ 2K [2 - 3 sin (X - ck'^)] 


(C-7) 


Q = 


W 

89 


1 H ^ (2-3 sin^ O' )1 sin^ O', 

. r\ ^ ^ Xn J i 

2|j,® 


+ 2K sin (X - Q'^) 


(C-8) 


S = 


3@^ 

89 


2 -1 

1 H ^ (2 -3 sin O') 

o -i 

2jji0 


- 2K sin 2(X - q-^) 


sin 2 O', 


(C-9) 


Steady-State Regression Rate 

The summation of the constant terms in Eq. (C-1) gives 


lif = - P cos cy 
o o 


(C-IO) 


and since P is inherently positive, f is negative for direct orbits 
(|o I < 90°). 
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Oscillatory Amplitudes and Frequencies 

Combination Terms . Those oscillations resulting from all three 
of the perturbing influences have the following amplitudes and fre- 
quencies . 


(JU. = \^ 

1 ^o 


Ai = - S cos O' 

1 o 


B, = - S cos lot 


o 


uo^ = 2\t 

2 ^o 


A„ = Q sin O' 
z o 


= 


A„ cos O' 

2 o 


o)^ = 20 - 2\if 

3 o 


A ^ . 2 o . 

A„ = - Q srn -r- sin o' 

3 2 o 


B^ = - A, 


(ju, = 2 0 - ^ 

4 o ^o 


2 ^ 

A, = S sin (1 + 2 cos cy ) 

4 2 o 


B, = - A, 


o) = 20 
5 o 


Ac = - P sin O' cos a 
5 o o 


B. - - A, 


CD , = 2 0 + ijr 


O ' o 

2 ^o 

A. = S cos -TT- (1-2 cos a ) 

D Z O 


B. = “ A, 


CD-, = 20 + 2ijr 

7 o 


o ^ 

2 o 


A-, = Q cos -T- sin cy 
7 2 o 


B-, = - A. 


Solar Terms . Those terms arising from the solar effect only are 

• 2 • 

as follows, where N represents ©^/Gq. 
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3 4 ^1 . 

- 7 " N cos -7T Sin O' 

4 z o 


3 4 

— N cos -rr sin la 
8 2 o 


3 2 

- -r N cos -rr sin o', cos ck' 

4 2 1 o 

0 !' 

3 2 1. 

y N cos — sin Of. cos 2o! 

4 2 1 o 


^10 


9 2 

- yr N sin o', sin 2a- 
16 1 o 


A 3 . 2 "1 . 

A, ^ - y N sin -y sin o', cos q- 

11 4 2 1 o 

3 2 ^1 

B, = - 7 “ N sin y- sin Oi. cos 2a 

11 4 2 1 o 


3 4 ^1 

At „= 7 N sin yy sin a 
12 4 2 o 

3 4 ""l 

y N sin y- sin 2ty 
12 o 2 o 


- 2(§ - 2i|; 


A 1 .2 o 

At - y N sin -y sin yr- sin a 
13 4 2 2 o 


» 13 " 
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=29 
14 o 


14 


- 2 «) - 

3 2 ^ X 2 

= 2^ N sin — sin cy^(l + cos 20^^) sin ~ 


14 


14 


R " 20 

15 o 


- 2 © 


15 


9 2 

— N sin O', sin 2a 
32 1 o 


Bi5 = - A 


15 


UD^ , = 29 
16 o 


16 


- 2 © + ij; 


3 2 ^1 

- — N cos ~ sin aj^(l 


9 

o \ 2 

2 cos <y ) cos — 
o z 


®16 = - A 


16 


0 ) 


17 


= 29 - 2@ + 2ijf 


17 


^ /O', ry a 

3 4 1 2 o . 

^ N cos cos -j* sin a 


Bi7 = - A 


17 


0 ) 


18 


20 + 2@ - 2T!f 

o ^ o 


O Of 

2 o 


3 4 '^l _ _ 

A- o " - 7- N cos sin -r- sin a 
18 4 2 2 o 


B 


18 


- A 


18 


(JO. 


19 


29 + 20 - \lf 

o ^o 


3 2 ^1 - w 

A^^ = - ^ N cos — sin a^(l + 2 cos a^) sin — 


o 0 / 

2 o 


®19 ” ■ *19 


“20 = + 24 

9 2 
A20 “ ■ J2 ^ 

®20 ' ■ ^20 
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are 


(JO, 


21 


20 + 2(9 + ijr 

o 

3 2 

A21 4 N sin -y sin cy^Cl 

®21 " ■ ^21 


o <y 

o . 2 o 

2 cos cy ) cos -rr- 
o 2 


uo, 


22 


20 + 20 + 2\l; 

o ^ o 

o /.O', ^ Of 

A _3 ,41 2o, 

A - 7 - N sin ~ cos — sin o- 
22 4 2 2 o 

R = - A 
22 ^22 


Lunar Terms , Those terms resulting entirely from the lunar effe 


* 2 

as follows, where M represents 0 /|J<6 : 

m o 


(JO, 


23 


^ - 2^ 
m ^o 


Aoo ~ TT M sin 2cy sin 2cy,, sin cy 
2J J2 m 1 o 


“ “ To M sin 2cy sin ry, sin 2cy 
2 j J2 ml o 


(JO, 


24 


t "I' 


^o/, “ ■ TT ^ 2cy [(cos ty + cos 2 q', ) cos cy 
24 io m^ 1 1 o 

- sin O'! sin cy ] 

1 o 


= TT ^ 20 - (cos O' + cos 2cy,, ) cos 2(y 

24 16 m l 1 o 


uo. 


25 “ ^m 
A 


25 


= 0 


B 


25 


■TT- M sin 20^ sin 2o', sin la 
32 m 1 o 
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OD, 


26 


t + t 

m 


^26 ” ' 2(y_[(cos <y^ - cos 2a^) cos 


16 


m 


sin O'- sin or 1 
1 o 


= - TT M sin 2<y (cos - cos 2or, ) cos 2a 
2d io mi i o 


0), 


27 


ik + 2\ir 

o 


^27 


- “ M sin 2 cy sin 2 q', sin o' 
32 m 1 o 


27 32 


M sin 2 q^ sin a-, sin 2 Qf 


m 


^o« = - 2i 


"28 

A, 


28 


B 


28 


4 ^1 


3 2 . * . 

- 77 M sin O' cos -r- sin o' 
o m 2 o 

3 2 4 ^1 

rrr M sin O' cos -^r^in 2a 
16 m 2 o 


^OQ - ilf 


"29 

A, 


29 


B 


29 


3 . 2 2 "^1 . 

- — M sin a cos “:r- sin or, cos a 

8 m 2 i o 

3 2 2 *^1 

— M sin a cos -tt sin a. cos 2a 

8 m 2 1 o 


^30 == 


^30 ^ 


'30 


9 2 . 2 . ^ 
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Appendix D 


OSCILLATORY DEVIATIONS FROM STEADY -STATE REGRESSION 
STATEMENT OF THE PROBLEM 

In the body of this Report it is shown that under the influ- 
ence of the assumed perturbing forces, a satellite orbit is subject to 
a steady-state regression of the orbital plane. This motion is repre- 
sented in Fig. 6 where the normal to the orbital plane traces out a 
circle on the reference sphere at a constant angular rate. 

Analytically, this motion is described by Eqs. (51) and (52), where 

it is seen that there are oscillatory terms superposed on the steady- 

state solutions for a and \|; . It is pointed out in Appendix C that these 

solutions for ca and tjf are valid only if the oscillatory terms are small. 

Therefore, this appendix investigates the magnitude of these residual 
oscillations . 

REFERENCE SYSTEM 

To describe this oscillatory effect it is convenient to define an 
X, Y coordinate system as shown in Fig. 26. This sytem moves in such 
a way that the XY plane remains perpendicular to the steady-state or- 
bital normal and also tangent to the reference sphere of Fig. 6. The 
origin thus lies at the point of tangency, and the X axis remains tangent 
to the steady-state Z axis trace of Fig. 6. In such a coordinate sys- 
tem, any oscillatory components in either o' or \|f results in a departure 
of the instantaneous orbital normal from its steady-state position at 
the origin of the X, Y system. 





Fig. 26 — Reference system for oscillatory deviations of the orbit normal 
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OSCILLATORY MOTION 


In the expressions for o' and \|r given by Eqs. (51) and (52), it is 
assumed that at zero time all of the oscillatory components are in phase 
This can be generalized by introducing a new time origin at a time t 
and measuring time from this. Thus, Eqs. (51) and (52) are modified as 
follows 
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By making t large, the phases of the various components are essentially 
random since the frequencies are assumed to be incommensurable. 

From Eqs. (D-1) and (D-2), the X, Y coordinates of the intersection 
of the instaneous orbital normal with the XY plane are given by 
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An examination of the expressions for the (ju^ values in Appendix C shows 
that some of these can vanish for particular combinations of inclination 
angle and orbital radius. Those frequencies which can become zero are 
as follows 
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The contours in the r O’ plane for which these frequencies vanish are 

o o 

shown in Fig. 27. 

When a particular does become very small, the amplitude of 
the corresponding low-frequency oscillatory terms in Eqs. (D-1) and 
(D-2) become large due to the 1/u)^ factor. However, the importance 
of such terms depends on the period of time over which the motion 
is observed. For this reason, in the determination of the effective 
oscillatory amplitude an observation period, T, is defined, and the 
magnitude of the effective oscillation depends on this value of T. 




Fig. 27 — Zero frequency loci 
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LEAS T-SQUARES FIT 

After selecting a value for the interval T, a linear least-squares 
fit, X, is determined for the function X over this interval in the form 

127 127 
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The amplitudes of sin cd^t and cos cd^t in Eqs . (D-6) and (D-7) are shown 
in Figs . 28 and 29 . 

The function X represents the average trend of the function X over 

T T 

the interval between ~ ^ ~ 2 ' mean-square deviation of X 

from X is then given by 
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Fig. 29 — Dependence of aj-j and b 
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In a similar manner the trend of the function Y can be determined as 
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where 
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As before, the mean-square deviation of Y from Y is given by 
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COMBINED EFFECT 

As a result of the X and Y variations described above, the total 
steady-state displacement of the normal of the orbital plane from the 
origin of the X, Y system is given by 



(D-18) 


Since this quantity is a constant for the interval T, it can be absorbed 

into the initial value of <y and the initial value of the regression 

o ° 

angle This amounts to shifting the origin of the X, Y system by 

an amount D . 

o 

However, there is still a steady drift rate of the orbital normal 
relative to this new origin, which can be expressed as 



(D-19) 


from which the total drift during the time interval T is given by 

AD = Dj^T (D-20)_ 

The mean-square deviation from this steady drift can be expressed 
in terms of the components as 

2 2 2 

*^D ^X "^Y 
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From Eqs . (D-9) and (D-17) the argument of Eq. (D-21) can be expressed 
as 
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Thus, the deviation of the orbital regression from a steady-state 
motion can be described by the drift, AD, and the deviation from this 
drift, 

NUMERICAL RESULTS 
Drift 

In the determination of AD it is seen from Eqs. (D-7) , (D-15), (D-19) 
and (D-20) that AD is a function not only of the averaging interval T but 
also the arbitrary time T. However, the result should be independent of 
T, and it appears reasonable to replace cos and sin by their root- 

mean-square value of 1//T in the evaluation of AD. 

By means of the relations developed above, the value of AD can be 
computed as a function of the orbital inclination angle of a synchronous 
altitude satellite for values of T equal to 50, 100 and 500 years. The 
results of these computations are shown in Fig. 30. It should be note*d 
that the summations are taken over only the twelve lowest frequency terms, 
which are ^23 *^32* found that the contributions 

of the terms whose frequencies involve 0^, @ and 0^, are negligible. In 
Fig. 30 it is seen that for values of T which are appreciably larger than 
the oscillatory periods included, the values for AD are independent of 



Fig.3G Drift angle as a function of Q!q (synchronous orbit) 


AD (deg) 
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T and less than .5°. However, as o' approaches 90° the value of or 

o 2 

f 

2\|;^ approaches zero. In this region of high inclination the summations 
in Eq. (D-19) are dominated by the terms and which have the form 

shown in Fig. 29. Thus, the fluctuations in Fig. 30 are simply a re- 
flection of Fig. 29 plotted as a function of ex . The magnitude of the 

o 

maximum at cx^ = 90 increases with the averaging time, T. 

For comparison, Fig. 31 shows the variation of AD for an orbital 

radius of 4 earth radii and for T equal to 100 years. In this case, 

several maxima exist, as would be expected from Fig. 27. The maximum 

at 7 50° results from both n) 2 ^ and approaching zero while 0 ) 2 ^ 

gives a maximum at o’ = 70°; as before, (D causes a maximum at = 90°. 

O’ z o 

The maxima again are due to the dominance of the terms associated with 
the near-zero frequencies , and in the case of those at 50° and 70° the 
form of Fig. 29 appears above and below the zero frequency position, re- 
sulting in a symmetrical peak at each of these positions. As in the 
previous case, these maxima are accentuated as T increases, while the 
general level of the rest of the curve is relatively constant for values 
of T greater than the oscillatory periods. 

On the basis of Figs. 30 and 31 as well as other cases not shown 
here, it is found that the general level of AD decreases as the orbital 
altitude decreases . ^ 

Mean-Square Deviation 

2 

In the determination of by means of Eqs. (D-21) and (D-22) it 

2 

is again found that depends on t as well as T. However, examination 
of the argument of Eq. (D-21) as expressed in Eq. (D-22) shows that the 
maximum value of X + Y is given by 



Fig .31 — Drift angle as a function of - 4 earth radii ) 
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Thus , 


an upper bound on can be established by the relation 

127 127 

< Y y (X + Y ) 
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The mean-square deviation is evaluated as a function of the orbital 
inclination of a synchronous altitude orbit for values of T equal to 50, 
100 and 500 years. The results are shown in Fig. 32. The summations of 
Eq. (D-24) are taken over the twelve lowest frequency terms, as in the 
case of AD. 

An examination of Fig. 32 shows that for all values of T the curves 

are essentially the same for the lower values of cy . However, as ry in- 

o o 

creases, goes through a maximum. The height of this maximum increases 

with T, while the corresponding value of approaches 90°. As in the 

case of AD this maximum results from the fact that ou^ is approaching zero 

and the summation in Eq. (D-24) is dominated by the term (X„o + Y^„) _ . 

22 22 max 

Figure 33 represents the limiting curve# for Eq. (D-22) when r = s. Thus, 

the heavy portions of these curves represent (X,, + Y.,) which appears 

11 ii^max 

2 

in Qjj. A comparison of Fig. 33 and Fig. 32 shows that the maximum in 

Fig. 32 is equivalent to a plot of Fig. 33 as a function of cy for i = 2. 

o 

In view of the long solution time required to determine the curves 
in Fig. 32, less detailed information is available for other orbital 
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Fig. 32 — Root-mean-square deviation of the orbit normal as a function of «, 

(synchronous orbit) 
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2 

altitudes. However, it is found that the magnitude of decreases 
with the orbital radius. 

2 

Figure 34 shows the details of the curve for an orbital radius 
of 4 earth radii and an averaging time of 500 years for inclination 
angles in the vicinity of cx^ = 50 . This corresponds to the region in 
which both and vanish. It is seen that maxima occur on either 

side of the value of for which the two frequencies vanish. The shape 
of these maxima is also related to that shown in Fig. 33 for i = 24 
and 27. 

DISCUSSION 

On the basis of the foregoing analysis, it is seen that the values 
of both AD and increase with orbital radius. Since the highest radius 
orbit considered here is that at synchronous altitude, Figs. 30 and 32 
represent upper bounds for AD and respectively. An examination of 
Figs. 30 and 32 shows that both the total drift and the root-mean-square 
deviation are less than .5 as long as the averaging time is appreciably 
larger than the periods of the oscillations included and if the orbital 
inclination angle, is less than about 75^. For orbital inclinations 

between 75° and 90°, significant peaks occur in the curves for both AD 
and . Thus, the steady-state motion described by Eqs. (57) and (58) 
is not valid at these high inclination angles where becomes very small. 
As indicated previously, the behavior of these high-inclination orbits 
is determined in Appendix F . 

Thus, it is seen that for inclination angles up to about 75° the 
normal to the orbital plane deviates less than .5° from its position 
as described by Eqs. (57) and (58). For an inclination greater than 
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Fig. 34 — Roof-mean-square deviation of the orbit normal as a function of a 
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75 , the regression may take place about the axis as described in 
Appendix F. Since these high-inclination orbits are of less interest, 
the contributions of the residual oscillatory terms for this second 
type of regression have not been evaluated. However, it does not appear 
that any large oscillatory contribution would result since many of the 
A, and B, values vanish at an inclination of 90°. 

L 1 




Appendix E 


LUNAR REGRESSION 


The satellite equations of motion developed in Appendix B should 
also apply to the case of the earth's most prominent satellite, the 
moon. However, as indicated earlier, it is necessary to go to a higher 
order approximation of the solution before the moon's orbital regression 
can be adequately described. The details of this determination are pre- 
sented below. 

MODIFICATION OF EQUATIONS OF MOTION 

The general equations of motion of an earth satellite, including per- 
turbations due to earth oblateness, the sun and the moon, are developed 
in Appendix B as Eqs . (B-57) through (B-60) . In applying these equations 
to the case of the moon, a number of simplif cations can be made. 

Of the three perturbations considered, the only significant one 
in the case of the moon is that due to the sun. At the moon's altitude 
the effect of earth oblateness is negligible, and the moon obviously 
cannot act as a perturbing influence on itself. Thus, the right side of 
Eqs. (B-57) through (B-60) should be modified by setting and 9^ equal 
to zero. 

In the original development, the reliction of the satellite on the ^ 
earth is neglected. However, this is not possible for a satellite 
as large as the moon. Thus, in Eq. (B-57) it is necessary to replace 
GMg by G(Mg + M^) . 

An examination of Eq. (54) shows that since the earth oblateness 
effect is negligible (J 2 = 0), the angle or^ is zei^ . Thus, the reference 
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plane becomes the plane of the ecliptic, and the x^, y^, and x^, y^, 
reference systems are identical. 

Finally, it can also be assumed that the moon's orbital inclination 
to the ecliptic is sufficiently small that cos cm is unity while sin ex 
is equal to a. 

Under the above assumptions, Eqs. (B-57) through (B-60) become 


^ _ 

dt^ 


puu = - 


(M^ 


+ M ) 
m 



(E-1) 


= 3p^@^(rj^ • i) (rj^ • j) (E-2) 

3@^(r^ • i) (r^ • k) sin 0 

oulf = — (E-3) 

' <y 

= 3©^(r^ * i) (r^ • k) cos 9 (E-4) 

where 

0) = 9 + \|f (E-5) 

and r has been replaced by p. 

By means of the direction cosines listed in Appendix A, the scalar 
products can be expressed as 

(r^ • i) = cos(0 - ® + i|;) 


(E-6) 
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• j) = 

- sin(9 - 0 + \|f) 

(E-7) 

(ri 

• k) = 

- O' sin(0 - i|;) 

(E-8) 


Substitution of Eqs. (E-6) through (E-8) in Eqs . (E-1) through (E-4) gives 
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(E-9) 


d , 2 , 
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(E-IO) 


(E-ll) 


a =- 


39^ 

uo 


>a cos 0 sin(0 - ij;) cos (0 + ifr - 0) (E-12) 


These expressions represent the desired equations of motion of the moon 
relative to the earth as affected by solar perturbations. 

METHOD OF SOLUTION 

In the case of the moon, it is not possible to determine ijf and d 
by substituting the unperturbed values of 0 and o) in Eqs. (E-ll) and 
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(E-12). Instead, it is necessary to take into account the effect of 
perturbations in ou and 0 on the equation for \|f and Q^, which in turn 
necessitates a solution for the perturbed in-plane motion as described 
by Eqs . (E-9) and (E-10). 


In-Plane Equations 

If CA)^ represents the observed sidereal rate of the moon around the 
earth, then the corresponding orbital radius in the absence of pertur- 
bations is defined by the relation 


P 

o 



(E-13) 


Substitution of Eq. (E-13) in Eq. (E-9) gives 


2 . ^ ^ • 2 

= . . . 9 . ° 1^1 + 3 cos(9 + i|f + 0)J 

dt p 


(E-14) 


as the form of the radial in-plane equation. 

While and uj^ are the unperturbed solutions of Eqs. (E-10) and 
(E-14). the quantities 6p and 6(A) represent the perturbations to these 
solutions which occur due to the sun. By making the following substi- 
tutions in Eqs. (E-10) and (E-14), 


p = p^ + 6p 

(E-15) 

U) = U) + 6(A) 

o. 

(E-16) 
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(JU 
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(E-17) 

(JL) t 
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(E-18) 
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the following perturbation equations are obtained in terms of 6p and 

6o). 
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The solution to these equations is given by 
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and 
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where m is the ratio of the earth's angular rate around the sun to the 

moon's angular rate around the earth, and 6p and 6u) are the initial 

o o 

values of 6p and 6 cjo. 

Since 6p^ and 6uu^ are arbitrary constants, they can be used to sim- 
plify the above solution by setting the coefficient of cos (v^t and the 
constant term in Eq. (E-23) equal to zero as follows 


26 p + 
o 


p6uo pm(5 + 4m) 

o o _ o 

(JD 12(1 - m) 


(E-24) 


36 p + 
o 


2p 6uu 
o o 

(JU 


2 2 
2p^m (3 - m - m ) 

3 - 8m + 4m^ 


(E-23^ 


These two conditions reduce Eqs. (E-22) and (E-23) to 


6p = - 


P m 


3p^m (2 - m) 


2(1 -m) (3 - 8m 4- 4m ) 


^ cos 2^t 
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and 

3ud^( 11 - 12m + 4m^) 

5uJ = ^ cos 2§t (E-27) 

4(1 - m) (3 - 8m + 4m ) 

Substitution of these values of 6p and 6ui) in Eqs. (E-15) and (E-16) 
gives the perturbed solutions for p and OJ as 

2 

P " ^o[^ ■ T ■ ^ m) cos 2itJ (E-28) 

^ ^ ^ iJ (E-29) 

where powers of m greater than three are neglected in the expansion of 
6 p and 6 u). 

From Eqs. (E-28) and (E-29) it is seen that as a result of the 

choice of 6 p^ and 6u)^, the perturbed solution still has the observed mean 

orbital rate . However, to achieve this, the mean orbital radius must 

2 

be less than the unperturbed radius by an amount p^m /6, which is of the 
order of 220 mi. 

Out-of-Plane Equations 

Regression . It is now possible to determine the moon's regression 
rate, i]; , by means of Eq. (E-11), which can be transformed into the follow- 
ing form: 

■ 2 

30 r 

ijf = - 1^1 - cos 20 - cos 2(0 - \jf) 


+ cos 2 ( 9 + ij( - 0) 


(E-30) 
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If y]i and 6i|f represent the steady-state and oscillatory components, re- 
spectively, of the regression rate, then can be represented as 

^ + 6ilf (E-31) 

and the problem reduces to a determination of i|f^ from Eq. (E-30). Be- 
fore this can be done, it is necessary to obtain expressions for the 
quantities 0 and ijf as functions of time as follows. 

From Eq. (E-5), 0 can be expressed as 

0 = uj - \lr 

= 0) - \lf + 6cjd - 6\!; (E-32) 

o ^o 

which can be integrated to give 

0 = (oj^ - il;^) t J 'S^dt - J 6i|fdt (E-33) 

Since the quantities 6uu and 6\|; by definition have no steady-state value, 
their integrals will represent small oscillatory angles. 

Similarly, the angle can be expressed as 

r • 

ijf = i]; t + I 6i|fdt (E-34) 

O V 

If Eqs. (E-17), (E-29), (E-32) and (E-34) are substituted in Eq. (E-30), 

4 

the following expression for is obtained to the order of m 
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The values of and 6\j; can be determined by successive iterations of 
Eq. (E- 35 ) as follows. As a first approximation 

. 2 

3 (ju m 119 
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while Is given by the remainder of Eq. (E- 35 ), neglecting the 
integral term 
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Integration of Eq. (E-37) gives 
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Substitution of Eq. (E-38) in Eq. (E-35) gives an improved determination 
of i|; in the form 
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3(jd m 
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1 II 2, 
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A second iteration is not necessary since no additional terms of the order 
4 

of m arise in the steady-state regression rate. Thus, the regression 
rate of the moon is given by 


= - 0 ) 
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(E-40) 


which is identical with the expression determined in Ref. 3 using the 
method of Delaunay. 

By equating the steady-state parts of Eq. (E -5) , the mean angular 

rate 9 is obtained as 
o 


(E-41) 
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which when combined with Eq. (E-40) gives 


9 
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UL) Fl 
OL 
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9 3 273 41 

32 ™ ■ 128 _ 


(E-42) 


Since the angle 0 is measured from the line of nodes, the period associ- 
ated with 9^ is the nodical period or the period between passages of 
the ascending node, while the period associated with co^ is the sidereal 
period as indicated previously. The relation of these two periods is 
thus obtained from Eq. (E-42) as 
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S 
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9 3 273 41 

32 m - i 2 S ^ _ 


(E-43) 


Inclination . It is also possible to determine the variation in 
orbital inclination by expanding Eq. (E-12) as follows 

. 2 

« 3@ (y r 

QC = - I sin 20 - sin 2(0 - @ + i|;) 

+ sin 2(® - ijf)J (E-44) 

A comparison with Eq. (E-30) shows that the oscillatory amplitudes in 
Eq. (E-44) are reduced considerably by the factor a which is equal to 
about .1 rad. Thus, the determination of <y is somewhat simpler than 
that of \[r described above. 

Equation (E-44) can be written in the form 
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+ sin 2(4 - tj (E-45) 


where o) is replaced by and the terms in Scu and S’]/ are neglected 
Integration of Eq. (E-45) gives 
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and since the terms on the right are small, Eq- (E-46) can be expressed 
to the order of m^ as 
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- cos 2((jo 


♦ o) 


+ cos 2(uo 


- 9) t] 


(E -47) 



153 



The first and largest of the oscillatory terms in Eq. (E-48) has an 

amplitude of 3m0' /8 or about 8.6 min of arc, while the other two have 
o 

amplitudes of about .6 min of arc. 

Thus, it is seen that the inclination of the moon's orbit remains 
essentially constant. 




155 


I 11 'a k=3%5^ fi B ^ ^3 



Appendix F 

REGRESSION OF HIGH INCLINATION ORBITS 

INTRODUCTION 

In the body of this Report, the solutions for the orbital inclination 
angle, cy, and orbital regression angle, tjf , are given by Eqs . (51) and (52). 
In order to obtain these expressions it was assumed that where ijr appears 
on the right side of Eqs. (37) and (38) it can be replaced by In 

view of the fact that is proportional to the cosine of the inclination 
angle, it becomes very small as the inclination approaches 90^. An exam- 
ination of the resulting solution for \jr given by Eq. (52) shows that under 
these conditions the oscillatory terms are no longer negligible. Thus, 
for high-inclination orbits, the above assumption does not hold. This is 
indicated by the increase in the values of AD and in the vicinity of 
9D° as shown in Figs. 30 and 32. While the assumption does appear to be 
valid up to inclinations of the order of 75 , it is of interest to in- 
vestigate the behavior of higher inclination orbits. 

METHOD OF ANALYSIS 

If the substitution of ^^t for \|f 
can be written in the form 

or = A^ sin \^ + A^ sin 2\^ 

^ H — IB, cos \1; 

^ ^o sin O' L 1 


is not made, then Eqs. (49) and (50) 
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+ / A, sin 0 ). t 

L t 1 


(F-1) 


i“3 


+ cos 2^ 
111 


^ B^ cos 


i*3 


+ 


(F-2) 



156 


While \|/ is also present in the other oscillatory terms, it constitutes 
a slowly varying phase angle of a higher frequency oscillation and can 
be neglected. However, when it stands alone, as in the first two oscil- 
latory terras, it must be considered. An examination of the numerical 
values of A and B shows that the summation terms in Eqs. (F-1) and (F-2) 
are negligible, so that 

O’ = A^ sin \lf + A^ sin 2\]/ (F-3) 


ijf = \lf + — 7 ^ 

^o Sin O' 


cos T|f + B^ 


cos 




(F-4) 


If the reference plane is defined in the same manner as before to make 
A^ and B^^ zero, then Eqs. (F-3) and (F-4) can be expressed as 

O' = Q sin O' sin 2i|f (F-5) 

ij; = - (P - Q cos 2\|r) cos a (F-6) 

where P and Q are defined by Eqs. (C-7) and (C-8) of Appendix C. 
Combination of Eqs. (F-5) and (F-6) gives the relation 


dcy _ Q sin oi sin 2\|f 
d\|f (P - Q cos 2i|f)cos ot 


which can be integrated to give 


(F-7) 


sin O' 


sin 01 

o 



Q cos 
Q cos 


2 | 

^ o 

2 \^ 


(F-8) 
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as the functional relation between O' and \|j . The quantities and 
represent the initial values of cy and , and, in accordance with Eqs . (51) 
and (52), is equal to zero when cn equals Thus, Eq. (F-8) becomes 


sin O' 


sin cy 

o 



---Q 

- Q cos 2ij; 


(F-9) 


However, for the purposes of this analysis it is more convenient to select 
the initial position when \{f is 90 and a has a value of Under these 

conditions, Eq. (F-8) becomes 


sin (X = sin QC 


1 

P + Q 

/p • 

- Q cos 2^ 


(F-10) 


while Q! and (X are related by the expression 
o o 


V P - Q 

FTq 


(F-11) 


By means of Eq. (F-10), the trace of the normal to the orbital plane on 
a unit sphere can be determined for any given value of In Fig. 35, 

several of these traces are shown for various orbital inclinations of a 
synchronous altitude orbit. In this figure, the angle ijf is measured in 
the 5 c^y^ plane from the negative y^ axis while cx is the elevation angle 
measured from the axis. It is seen that for values of (x^ up to about 
75*^ the traces encircle the z^ axis at a relatively constant value for 
the inclination angle, ry . This corresponds to the regression described 


in the body of the report and pictured in Fig. 6. 
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Normal to reference plane 
(stable axis) 



Fig. 35 — Regression of a synchronous altitude orbit 
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For values of a* such that 
o 

the nature of the traces change and they become elongated ellipses which 
encircle the axis. 

The boundaries between the two types of trace correspond to the con- 
ditions 

or a value = 78°53' for synchronous altitude orbits. The actual con- 
tours lie in planes including the axis and inclined at angles of 
±11 7 to the x^y^ reference plane. 

The periodicity of the motion corresponding to these traces can be 
determined if the angle of is eliminated between Eqs . (F-6) and (F-10) to 
give 

t (P - Q cos 2i|f) (P cos - Q sin “ Q 2^) (F-14) 


which can be solved for t in the form 


T 

= - J 

n 


dx 


- Q COS 2\jf) (P cos^ - Q sin^ - Q cos 2ijf) 


(F-15) 


By means of the transformations given in Ref. 4, Eq. (F-15) can be ex- 
pressed as an elliptic integral of the first kind. The resulting 
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expression for the period of the motion is given by the following equations. 
If -sin + Q ’ then 


T 



COS 


(y 

o 




d$ 


k Sin 


where 


k 


2 



2 

tan 


/ 



(F-16) 


(F-17) 


If sin (y > 




L . - . 

P + Q 


T = 


^J^Q(P + Q) sin a 


t 


d$ 




,2 .2 

k sin 


(F-18) 


where 


k 


2 


? - Q 

20 tan ot' 
^ o 


(F-19) 


Thus, it is seen that the period given by Eq. (F-16) corresponds to the” 
more conventional regression about the axis while that given by Eq. (F-18) 
is associated with the elongated elliptical traces of Fig. 35. A plot of 
the regression period as a function of is shown in Fig. 36. 
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Orbital inclination oto (deg) 

^ = 90 ° 


Fig. 36 — Periods of regression about X| and zj axes 
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RESULTS AND DISCUSSION 

In the body of the Report it is shown that the plane has the 

property that an orbit in this plane would maintain its orientation 
relative to inertial space. However, an examination of Fig. 35 shows 
that a second stable configuration exists with the orbit in the 
plane. Since the earth's axis lies in this plane, such an orbit would 
not only remain fixed relative to inertial space but would remain polar 
relative to the earth. On the other hand, an orbit established in the 
^1^1 unstable equilibrium and may regress about either the 

Xj^ or axis, depending on the direction of its initial disturbance. 

In order to compare the regression periods determined by this method 
with those shown in Fig. 7, it is necessary to express Eq. (F-16) in terms 
of the inclination angle when ^ is equal to zero. This is done by com- 
bining Eqs . (F-11) and (F-16) to give 




determined in the body of the Report as Eq. (62) . Figure 37 is a plot 
of Eqs. (F-20) and (F-21) for synchronous altitude orbits. It is seen 
that the agreement is excellent up to large values of at . 
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{^ = 0 ) 


Fig. 37 Comparison of Eqs, ( F - 20 ) and ( F - 21 ) for regression period 
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The maximum variation of orbital inclination during regression about 
the z, axis is equal to the difference between or and ty^ . This can be 
seen in Fig. 35, where the values are indicated relative to the actual 
position of the trace when it crosses the y-^-^ plane. It is seen that 
the variation is relatively small, having reached a value of only 5° at 
an inclination of 75°. Thus, the representation of this type of re- 
gression by a conical motion of the normal to the orbital plane as shown 
in Fig. 6 also appears to be valid as long as the regression is about 
the axis. 

The results obtained in this appendix are in excellent agreement 
with those obtained in Ref. 5 by Allan and Cook who used eigenvalue 
methods to determine the nature of the regression and its periodicity. 
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